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Foreword 


The  Government  of  Canada  and  the  provincial  and  territorial  gov- 
ernments are  working  with  the  agriculture  and  agri-food  industry 
and  interested  Canadians  to  develop  an  architecture  for  agricul- 
tural policy  for  the  21st  century.  The  objective  of  the  Agricultural 
Policy  Framework  (APF)  is  for  Canada  to  become  the  world  leader 
in  food  safety  and  food  quality,  innovation  and  environmen- 
tally-responsible production.  To  contribute  to  these  goals,  Agricul- 
ture and  Agri-Food  Canada  (AAFC)  has  an  ongoing  research 
program  to  provide  information  on  the  effects  of  agricultural  pol- 
icy and  technology  scenarios  on  the  environment  and  on  the  eco- 
nomic performance  of  the  agriculture  sector. 

Included  in  this  work  program  is  a  project  to  improve  our  farm 
level  data  on  the  cost  of  production  and  the  farm  management 
practices  for  economic  and  environmental  analysis.  As  part  of  this 
effort  to  improve  our  data,  this  report  evaluates  an  analytical 
method,  called  Maximum  Entropy  (ME),  for  its  effectiveness  in 
extracting  detailed,  enterprise-level  cost  of  production  information 
from  whole-farm  data  sets.  The  ME  method  has  been  shown  to  be 
a  promising  and  cost-effective  option  for  obtaining  these  enter- 
prise-level estimates  from  whole  farm  data  sets  already  available. 
Traditionally,  direct  collection  of  these  estimates  has  been  difficult 
and  costly.  Further  research  on  the  ME  methodology  is  underway. 
This  study  will  improve  our  capacity  to  carry  out  analysis  in  sup- 
port of  the  APF  by  expanding  our  understanding  of  the  linkages 
between  farm  operations  and  environmental  indicators. 

Any  policy  views,  whether  explicitly  stated,  inferred  or  interpreted 
from  the  contents  of  this  report  do  not  necessarily  reflect  the  views 
or  policies  of  AAFC. 
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Executive  summary 


Governments  have  gone  to  considerable  lengths  in  the  past  to 
obtain  economic  data  on  the  individual  enterprises  of  multi-output 
farms.  For  policy  purposes,  governments  need  information  with 
respect  to  the  costs  per  unit  of  output,  measured  by  so-called 
input-output  coefficients,  or  (preferably)  cost-allocation  coeffi- 
cients. However,  collecting  this  information  is  costly,  as  considera- 
ble resources  have  been  devoted  to  apportion  the  whole-farm 
accounting  data  among  the  individual  enterprises.  Consequently, 
there  is  a  need  for  policy  analysts  and  practitioners  to  access  a 
methodology  that  can  offer  reliable  estimates  at  a  significantly 
lower  cost. 

The  Maximum  Entropy  (ME)  methodology  is  a  valid  option  for 
obtaining  these  enterprise-level  estimates  from  whole-farm  data. 
The  ME  estimates  have  the  desirable  property  of  being  consistent 
with  the  available  data  and  known  constraints  while,  at  the  same 
time,  not  requiring  a  host  of  assumptions  by  the  analyst.  As  well, 
since  the  estimates  of  the  underlying  commodity  cost  (and  hence, 
production)  structures  are  obtained  from  data  sets  already  availa- 
ble, such  as  the  Taxfiler  Database  or  the  Net  Income  Stabilization 
Account  (NISA)  Database,  the  results  can  be  used  to  improve  the 
Canadian  Regional  Agricultural  Model  (CRAM)  and  other  farm- 
level  models  without  undertaking  more  costly  data  collection. 

The  purpose  of  this  study  is  to  demonstrate  the  relevance  of  the 
ME  methodology  in  estimating  cost-allocation  coefficients  and  to 
implement  the  ME  methodology  using  a  data  set  from  a  sample  of 
Saskatchewan  crop  farms.  Furthermore,  the  study  assesses  the 
"quality"  of  the  ME  methodology  by  looking  at  the  precision  with 
which  the  enterprise-level  cost  structure  can  be  extracted  from  the 
whole-farm  data.  Specifically,  the  estimation  results  for  a  sample 
of  farms  are  compared  to  the  actual  enterprise-level  data  that  are 
available  for  the  same  sample  of  farms. 
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Executive  summary 


Two  applications  of  the  ME  methodology  are  conducted.  The  first 
application  estimates  cost  allocation  coefficients  assuming  that 
each  cost  item  is  a  linear  function  of  several  output-specific  reve- 
nues. The  ME  methodology  is  then  applied  to  a  complete  set  of  lin- 
ear relations,  linking  each  cost  item  to  several  output  variables.  In 
addition,  the  empirical  model  is  adjusted  to  take  into  account  sev- 
eral "zero  observations"  where,  for  some  farms,  the  sample  expen- 
ditures on  certain  input  categories,  including  net  operating 
income,  are  equal  to  zero  or  are  negative.  The  second  application 
includes  the  first  application  but  also  takes  into  account  the  partic- 
ular nature  of  Saskatchewan  agriculture,  in  which  a  significant 
proportion  of  land  is  put  aside  each  year  due  to  adverse  agro-cli 
matic  conditions.  Maintaining  this  so-called  "fallowed  land 
requires  certain  farm  operations,  including  tillage  and  weed  con- 
trol. Since  these  operations  have  significant  costs,  the  costs  need  to 
be  considered  in  the  analysis  by  adjusting  the  ME  methodology 
accordingly. 


it 


Empirical  implementation  of  the  ME  methodology  is  conducted 
using  the  General  Algebraic  Modelling  System  (GAMS)  computer 
program.  Cost-allocation  coefficients  are  estimated  for  five  crop 
enterprises:  wheat,  other  grains,  canola,  other  oilseeds  and  other 
crops.  From  the  empirical  results,  two  main  findings  emerge: 

The  ME  methodology  is  a  valid  approach  to  estimate  cost-alloca- 
tion coefficients,  in  the  sense  that  the  estimated  coefficients  fall 
within  the  range  of  their  observed  counterparts. 

Some  minor  qualifications  are  needed,  depending  on  the  category 
of  crops  under  consideration.  The  estimated  and  observed  cost- 
allocation  coefficients  are  quite  comparable  for  "Wheat,"  "Other 
Grains,"  and  to  a  lesser  extent  "Canola."  However,  for  the  two 
remaining  crop  categories,  "Other  Oilseeds"  and  "Other  Crops," 
some  caution  is  in  order  when  comparing  the  estimated  and 
observed  cost  allocations.  The  discrepancies  found  are  not  really 
surprising  given  the  heterogeneity  in  those  categories.  Neverthe- 
less, further  investigation  may  be  necessary  by,  for  example,  creat- 
ing more  homogenous  crop  categories. 

Given  the  exploratory  nature  of  the  present  study,  further  work  is 
still  needed  to  obtain  more  efficient  and  realistic  estimates-for 
example,  by  using  larger  samples,  exploiting  the  "richness"  of 
panel  data  and/ or  employing  more  reliable  prior  information  on 
the  support  values  of  the  cost-allocation  coefficients.  Some  other 
practical  and  methodological  problems  remain.  In  particular,  one 
issue  that  needs  further  investigation  is  how  to  model,  in  a  more 
satisfactory  way,  the  decision  process  of  fallow  land  maintenance. 
It  would  require  the  use  of  sample  data  for  several  years  to  capture 
the  delayed  effect  of  fallow  on  the  returns  of  crop  farmers,  or  the 
use  of  more  homogeneous  crop  farms-for  example,  crop  farms  that 
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are  located  in  the  same  region  or  that  are  using  similar  technolo- 
gies. Another  issue  that  deserves  more  emphasis  is  the  necessity 
that  crop  and  livestock  enterprises  be  described  and  defined  with 
greater  accuracy.  Finally,  it  is  necessary  to  raise  the  issue  of  hetero- 
scedasticity,  because  input  structures  are  likely  to  vary  according 
to  the  volume  of  production  on  each  farm.  This  issue  could  be 
examined  in  future  research  using  appropriate  model  specifica- 
tions. 
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Section  1:   Introduction 

In  conducting  an  economic  analysis  into  farm  level  response  to  changes  in  markets,  technol- 
ogy or  policy,  a  clear  understanding  of  the  production  processes  and  the  associated  cost 
structures  is  required.  Unfortunately,  this  type  of  detailed  data  is  rarely  available,  as  farmers 
do  not  usually  keep  this  type  of  management  or  financial  records  or  if  they  do,  it  is  not  avail- 
able to  those  conducting  the  economic  analysis.  Normally,  the  source  of  data  that  is  readily 
available  is  whole-farm  data,  in  which  the  revenue  obtained  from  various  enterprises  may  be 
specified,  but  costs  for  the  various  inputs  are  lumped  together  for  all  enterprises. 

Governments  have  gone  to  considerable  lengths  in  the  past  to  obtain  economic  data  on  the 
individual  enterprises  of  multi-output  farms.  For  policy  purposes,  governments  need  infor- 
mation with  respect  to  the  costs  per  unit  of  output,  measured  by  so-called  input-output  coeffi- 
cients or  (preferably)  cost-allocation  coefficients.  However,  collecting  this  information  is  very 
costly,  as  considerable  effort  has  to  be  expended  to  apportion  the  whole-farm  accounting 
data  among  the  individual  enterprises.  As  a  result,  there  is  a  need  for  policy  analysts  and 
practitioners  to  access  a  methodology  that  can  provide  reliable  estimates  at  a  significantly 
lower  cost. 

The  Maximum  Entropy  (ME)  methodology  is  a  valid  option  for  obtaining  these  enterprise- 
level  estimates  from  whole-farm  data.  The  ME  estimates  have  the  desirable  property  of  being 
consistent  with  the  available  data  and  known  constraints  while,  at  the  same  time,  not  requir- 
ing a  host  of  assumptions  by  the  analyst.  As  well,  since  the  estimates  of  the  underlying  com- 
modity cost  (and  hence,  production)  structures  are  obtained  from  data  sets  already  available, 
such  as  the  Taxfiler  Database  or  the  Net  Income  Stabilization  Account  (NISA)  Database,  the 
results  can  be  used  to  improve  the  Canadian  Regional  Agricultural  Model  (CRAM)  and  other 
farm-level  models  without  undertaking  more  costly  data  collection. 

The  objective  of  this  study  is  to  use  data  that  currently  exist  within  the  Research  and  Analysis 
Directorate  (RAD)  of  Agriculture  and  Agri-Food  Canada  (AAFC)  to  demonstrate  and  to  test 
the  proposed  ME  methodology.  Furthermore,  the  study  will  assess  the  "quality"  of  the  ME 
methodology  by  looking  at  the  precision  with  which  the  enterprise-level  cost  structure  can  be 
extracted  from  whole-farm  data.  Specifically,  the  estimation  results  for  a  sample  of  farms  will 
be  compared  to  the  actual  enterprise-level  data  that  are  available  for  the  same  sample  of 
farms. 
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This  study  describes  two  applications  of  the  ME  methodology  to  a  data  set  from  a  sample  of 
Saskatchewan  crop  farms.  The  first  is  a  direct  application  of  the  ME  methodology,  devel- 
oped by  Leon  et  al.  (1999)  for  a  group  of  farms  in  the  Brittany  Region  of  France,  where  each 
cost  item  is  assumed  to  be  a  linear  function  of  several  output-specific  revenues.  The  ME 
methodology  is  then  applied  to  a  complete  set  of  linear  relations,  linking  each  cost  item  to 
several  output  variables.  In  addition,  the  empirical  model  is  adjusted  to  take  into  account 
several  "zero  observations"  where  for  some  farms,  the  sample  expenditures  on  certain  input 
categories,  including  net  operating  income,  are  equal  to  zero  or  are  negative.  The  second 
application  includes  the  first  application  but  also  takes  into  account  the  particular  nature  of 
Saskatchewan  agriculture,  in  which  a  significant  proportion  of  land  is  put  aside  each  year 
due  to  adverse  agro-climatic  conditions.  Maintaining  this  so-called  "fallowed  land"  requires 
certain  farm  operations,  including  tillage  and  weed  control.  Since  these  operations  have  sig- 
nificant costs,  the  costs  need  to  be  considered  in  the  analysis  and  the  ME  methodology  has  to 
be  adjusted  accordingly. 

The  present  study  is  based  on  1994  accounting  data  from  a  sample  of  38  Saskatchewan  crop 
farms,  as  supplied  by  the  management  of  the  Farm-Level  Database  System  of  the  Strategic 
Policy  Branch,  AAFC.  The  data  are  part  of  the  Top  Management  Model2  developed  in  the 
early  1980s  by  Richard  Schoney  at  the  University  of  Saskatchewan  (Schoney  1991). 

All  the  farms  in  the  sample  produce  several  crops,  including  grains  (wheat,  barley,  oats  and 
rye),  oilseeds  (canola,  linseed  and  mustard  seed),  specialty  crops  (canary  seeds,  caraway 
seeds,  lentils,  etc.)  and  forage  crops.  From  a  total  of  38  farms  in  the  sample,  eight  farms  were 
excluded  from  the  analysis  because  they  were  involved  either  in  single-crop  production 
(wheat)  or  in  non-crop-related  activities,  such  as  seed  cleaning  and  custom  work.  In  other 
words,  only  multi-crop  farms  producing  grains,  oilseeds,  specialty  crops  and  forage  crops 
were  considered  in  the  analysis.  In  addition,  the  farms  in  the  sample  are  not  concentrated  in 
one  specific  region  of  Saskatchewan  but  are  rather  evenly  distributed,  thus  providing  a  bal- 
anced representation  of  its  agricultural  sector. 

The  remainder  of  this  study  is  organized  into  four  sections.  Section  2  provides  a  preliminary 
overview  of  the  Saskatchewan  farm.  Section  3  presents  the  ME  methodology.  Section  4 
reports  the  two  applications  and  we  discuss  the  corresponding  results.  Section  5  is  a  sum- 
mary of  the  main  findings,  along  with  recommendations  for  further  research.  The  list  of  ref- 
erences and  five  appendices  on  the  data  and  selected  crop  categories  appear  at  the  end  of  the 
study. 


The  initial  purpose  of  this  study  was  to  apply  the  ME  methodology  to  two  groups:  a  group  of  crop  farms  in 
Saskatchewan  and  a  group  of  livestock  farms  in  Ontario.  Eventually  however,  only  the  Saskatchewan  crop 
farms  were  retained,  because  the  sample  of  Ontario  livestock  farms  was  not  suitable  for  the  present 
analysis.  Most  of  the  farms  in  the  Ontario  sample  were  characterized  by  an  overly  strong  specialization 
toward  hogs  (no  multi-enterprises  farms). 

The  Top  Management  Model  was  initially  designed  as  a  farm  management  extension  tool,  which  provided 
highly  sophisticated  enterprise-level  and  farm-level  financial  analysis  of  the  current  and  projected  periods. 
Later,  the  model  was  adapted  to  serve  as  a  data-collection  tool.  The  Saskatchewan  data  were  collected 
through  the  Top  Management  Workshops,  a  series  of  extension  workshops  held  annually  from  1987  to 
1994.  The  workshops  were  conducted  by  the  University  of  Saskatchewan  and  received  joint  funding  from 
AAFC  and  the  Saskatchewan  Department  of  Agriculture  and  Food. 
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Section  2:   Preliminary  data  analysis 

A  preliminary  data  analysis  was  necessary  for  three  reasons.  It  enabled  us  to  organize  the 
data  according  to  specific  cost  and  revenue  categories,  to  correct  several  data  inconsistencies, 
and  to  ensure  the  accounting  balance  between  the  revenues  and  costs  for  each  farm. 

In  Section  2,  we  discuss  several  aspects  of  the  data  sample  of  Saskatchewan  farms.  We  define 
the  various  cost  and  revenue  items  retained  for  the  ME  estimation  of  cost-allocation  coeffi- 
cients in  Section  2.1.  Then  we  present  the  main  characteristics  of  the  selected  Saskatchewan 
crop  farms  in  Section  2.2,  followed  in  Section  2.3  by  a  discussion  of  agro-climatic  conditions 
prevailing  in  Saskatchewan.  We  review  the  inconsistencies  appearing  in  the  data  sample  of 
Saskatchewan  crop  farms  and  provide  ways  (adjustments)  adopted  to  overcome  them. 

2.1     Aggregate  cost  and  revenue  categories 

To  make  the  estimation  of  the  cost-allocation  coefficients  feasible,  the  original  farm  data  were 
aggregated  according  to  a  limited  number  of  categories  (see  Box  1).  Specifically,  five  aggre- 
gate revenue  (output)  and  nine  aggregate  cost  (input)  categories  were  considered.  The  selec- 
tion of  these  categories  was  carried  out  in  consultation  with  experts  from  AAFC.  The  selected 
categories  reflect  the  major  characteristics  of  Saskatchewan  farms  in  terms  of  crop  orientation 
and  production  systems. 

Each  of  the  aggregated  output  categories  comprises  a  wide  range  of  crop  enterprises,  which 
are  listed  in  Appendix  A. 

In  addition  to  the  nine  input  categories,  the  analysis  considers  "fallow  costs"  (see  Footnote  24 
for  details  on  the  definition  of  fallow  costs).  "Fallow  land"  was  considered  here  which 
included  not  only  fallow  land,  but  also  the  land  planted  with  crops  aimed  at  improving  soil 
fertility.  The  crops  of  concern  are  mainly  green  manure  and  green  feed.  Appendix  B  contains 
the  figures  on  net  revenues  and  costs  for  each  farm  included  in  the  sample. 
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Box  1:  Overview  of  revenue  and  cost  categories 

Revenues  (output) 

Costs  (input) 

Wheat-excluding  wheat  grown  for  seed  production3 

Application  1 

Other  Grains-including  (feed  and  malt)  barley,  oats  and  ryeb 

Seeds 

Canola 

Fertilizers 

Other  Oilseeds-including  flax,  sunola  and  mustard  seeds0 

Pesticides 

Other  cropsd 

Other  direct  material  inputs 

Fuel,  power  and  grease 

Repairs 

Paid  salaries 

Other  fixed  cash  expenses 

Net  operating  income 

Application  2  also  includes  fallow  costs 

a.  Spring  wheat  and  durum  wheat  grown  for  seed  production  had  to  be  excluded  from  the  analysis, 
because  these  two  types  require  special  treatment.  Wheat  grown  for  seed  production  corresponds  to 
crop  enterprises,  such  as  "certified  and  pedigree  wheat,"  which  are  viewed  as  producing  specialty 
crops  and  hence,  are  part  of  "Other  crops," 

b.  Although  viewed  as  a  hybrid  grain,  triticale  was  categorized  as  a  specialty  crop,  and  was  not  included 
in  "Other  grains."  This  crop  plays  a  minor  role  in  the  present  analysis,  as  it  is  grown  only  by  a  small 
number  of  farms. 

c.  Similar  to  wheat,  some  oilseeds  are  grown  for  seed  production.  In  this  case,  they  are  considered  a 
specialty  crop  and  hence,  are  included  in  "Other  crops." 

d.  This  category  includes  all  the  crops  excluded  from  the  other  categories.  It  is  made  up,  mainly,  of 
specialty  crops  and  forage  crops. 

2.2     Farm  characteristics 

The  "average"  farm  in  the  sample  earns  a  gross  operating  income  of  $258,000  on 
2,186  hectares  (see  Table  1).  Wheat  and  canola  are  the  two  major  crops  in  the  sample  repre- 
senting 33  percent  and  31  percent  of  total  revenues,  and  33  percent  and  19  percent  of  total 
hectares,  respectively.  The  remaining  crops  represent  about  one  third  of  gross  farm  income. 
The  net  operating  income  of  the  average  farm  is  about  $98,000,  whereas  total  costs  represent 
62  percent  of  total  gross  revenue.  A  closer  look  at  the  various  costs  reveals  that,  except  for 
other  fixed  cash  expenses  at  16  percent,  each  represents  less  than  10  percent  of  total  costs. 

The  farms  in  the  sample  show  considerable  heterogeneity  in  acreages,  revenues  and  costs  as 
shown  by  the  statistical  indicators  presented  in  the  last  three  columns  of  Table  1.  Specifically, 
the  coefficients  of  variation  are  quite  large  (larger  than  100  percent).  This  finding  is 
confirmed  by  the  magnitudes  of  the  range  (minimum-maximum)  of  the  revenue  and  cost 
categories.  However,  the  farms  show  considerable  variation  in  net  operating  income.  While 
some  farms  experienced  a  negative  net  operating  income,  the  most  profitable  farms  have  a 
net  operating  income  that  is  four  times  greater  than  the  average  value. 
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2.3  Summerfallow  and  crop  rotations 

The  agro-climatic  conditions  prevailing  in  Saskatchewan  require  summerfallow  and  specific 
crop  rotations.  The  crop  rotations  vary  according  to  three  soil  zones  that  run  through  the 
province:  southern  brown  soil  (SBS),  central  dark  brown  soil  (CDBS)  and  black  soil  (BS). 
Summerfallow  plays  an  integral  part  of  the  crop  rotation  system  in  the  SBS  zone  and  in  the 
drier  parts  of  the  CDBS  zone.  Farmers  in  the  BS  zone  and  parts  of  the  CDBS  zone,  on  the 
other  hand,  have  more  or  less  discontinued  the  use  of  summerfallow.  Given  that  the  BS  and 
CDBS  zones  receive  greater  precipitation  and  have  higher  soil  fertility,  the  adoption  of  new 
production  techniques,  tillage  equipment  and  crop  varieties  has  gradually  allowed  for  con- 
tinuous cropping.4  The  sample  data  clearly  reflect  these  characteristics  of  the  province's  agri- 
culture and  the  underlying  agro-climatic  conditions.  Six  farms  in  the  sample  do  not  use 
summerfallow  which,  for  the  other  farms,  represents  on  average  about  one  fifth  of  the  total 
acreage. 

2.4  Data  inconsistencies  and  appropriate  adjustments 

Preliminary  inspection  of  the  original  data  set  revealed  inconsistencies  that  required  appro- 
priate adjustments.  These  inconsistencies  stem  partly  from  the  fact  that  the  original  data 
were  recorded  according  to  two  different  methods.  Specifically,  revenues  (and  costs)  were 
recorded  on  either  a  cash  basis  or  an  accrual  basis.  The  two  methods  are  different  only  in  how 
cash  costs  and  sales  are  recorded:  the  former  distinguishes  cash  costs  from  farm  input  and 
crop  inventory  changes,  while  the  latter  records  costs  of  the  inputs  effectively  used  by  the 
farms  and  crop  outputs  produced  by  the  farms. 

Given  the  purpose  of  the  present  study  to  estimate  cost-allocation  coefficients  associated  with 
the  inputs  effectively  used  by  the  Saskatchewan  farmers,  the  accrual  cost  data  approach  was 
adopted.  This  choice  is  more  in  line  with  the  economic  notion  of  costs  (and  revenues).  Also, 
the  database  reporting  costs  and  revenues  on  an  accrual  basis  was  organized  for  each  crop 
enterprise,  which  allowed  for  assessing  the  performance  of  the  ME  methodology.  By  con- 
trast, the  data  file  presenting  costs  and  revenues  on  a  cash  basis  provides  only  totals,  with  no 
breakdown  by  categories  of  crops. 


Typical  rotations  in  these  regions  are  #1  Fallow-wheat-wheat,  #2  Fallow-wheat-feed  grain  (barley,  oats  and 
rye),  #3  Fallow-flaxseed  or  mustard  seed-wheat-feed  grain,  #4  Fallow-wheat-lentils-wheat,  and  #5  Fallow- 
wheat-feed  grain-forage  crop  (three-five  years).  The  length  of  the  rotation  depends  on  the  agro-climatic 
conditions  of  the  local  area.  Dry  regions  with  low  soil  fertility  have  short  rotations  (e.g.  #1  and  #2)  while 
areas  with  more  moisture  and  greater  fertility  have  longer  rotations  (e.g.  #3  and  #4). 

Typical  rotations  in  these  regions  include  #1  Wheat-wheat-feed  grain,  #2  Canola-peas-wheat-wheat, 
#3  Canola-feed  grain-wheat-peas-wheat,  #4  Wheat-feed  grain-forage  crop  (for  two-four  years). 
Summerfallow  may  still  be  used  periodically  to  treat  disease,  insect  and  weed  problems  or  during  drought 
years,  but  is  no  longer  a  regular  component  of  the  rotation. 
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Table  1:  Saskatchewan  farm  data 


Mean 

c.v. 

Min 

imum 

Maximum 

(ha) 

($) 

(percent) 

(percent) 

(ha) 

($) 

(ha)              ($) 

Acreages 

Fallow  land 

362 

166 

929 

0 

1,229 

Wheat 

723 

33.2 

1159 

0 

4,640 

Other  grains 

201 

93 

107.4 

0 

805 

Canola 

421 

18.9 

121  4 

0 

2,320 

Other 
oilseeds 

102 

47 

240  1 

0 

1,300 

Other  crops 

377 

17.3 

84.1 

0 

981 

Total 

2.186 

100.0 

82.0 

0 

10,040 

Revenues 

Wheat 

85.937 

33  3 

94.4 

0 

387,450 

Other  grains 

26,603 

10.3 

136  3 

0 

131,800 

Canola 

79.628 

30.9 

1257 

0 

372,938 

Other 
oilseeds 

11,001 

4.3 

2330 

0 

132,600 

Other  crops 

54,538 

21.2 

137.3 

0 

367,590 

Total 

257,707 

100.0 

899 

65,969 

1,025,226 

Costs 

Seeds 

14,727 

5.7 

867 

981 

59,645 

Fertilizers 

22,194 

86 

97.8 

0 

95,770 

Pesticides 

24,346 

9  5 

116.7 

4.325 

145,437 

Other  direct 

material 

11,584 

45 

70.2 

1,450 

32,309 

inputs 

Fuel,  power 
and  grease 

14,584 

5  7 

844 

0 

62,670 

Repairs 

14,491 

56 

86.2 

0 

66,200 

Paid 
salaries 

15,759 

6  1 

173.1 

0 

138,264 

Other  fixed 

cash 

42,034 

16.3 

82.3 

2,510 

172,018 

expenses 

Total 

159,719 

62.0 

828 

29,940 

670,550 

Net 

operating 

97,988 

380 

1168 

-12,031 

406,514 

Income 

Notes:  C.V.  =  Standard  deviation*100/Mean 

The  last  two  columns  present  the  minimum  and  maximum  of  each  item  obtained  from  the  data  sample 
of  Saskatchewan  farms  and  for  this  reason  all  revenue,  cost  and  net  operating  income  items  do  not  add 
up. 
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Before  applying  the  ME  methodology,  one  has  to  ensure  that  each  cost  item  be  measured 
similarly  by  using  the  accrual  and  cash  cost  methods.  To  do  so  required  the  reconciliation 
and  harmonization  of  the  two  databases  (supplying  costs  on  either  a  cash  or  an  accrual  basis, 
respectively).  Some  cost  categories,  which  were  recorded  for  some  farms  as  custom  work  in 
the  cash  cost  database,  were  recorded  differently  in  the  accrual  cost  database.  This  method  of 
recording  is  certainly  valid  for  expenses  associated  with  seed,  fertilizer  and  pesticide 
operations.  It  was  straightforward  to  reconcile  all  the  accrual  and  cash  cost  data  associated 
with  direct  material  inputs  (i.e.  seeds,  fertilizers,  pesticides,  other  direct  material  inputs); 
fuel,  power  and  grease;  and  repairs.  On  the  other  hand,  for  paid  salaries  and  other  fixed  cash 
expenses,  it  was  impossible  to  reconcile  the  accrual  and  cash  cost  data.  In  fact,  too  many 
discrepancies  existed,  and  there  was  no  way  that  a  thorough  investigation,  aimed  at  reconcil- 
ing cost  data  for  these  two  categories,  could  be  conducted.  For  this  reason,  these  two  cost 
categories  were  measured  by  using  their  cash  cost  proxies.  Finally,  the  net  operating  income 
was  obtained  residually  by  taking  the  difference  between  total  accrued  sales,  as  measured  by 
total  net  returns  in  the  crop  enterprises  database  and  the  total  cost. 

Another  inconsistency,  which  is  linked  to  crop  activities,  stems  from  the  fact  that  many  farms 
in  the  sample  are  also  involved  in  some  non-crop  activities,  such  as  seed  cleaning,  custom 
work,  and  non-agricultural  activities  (e.g.  bird  habitat  conservation).  It  is  difficult  to  consider 
these  farms,  because  there  is  no  way  to  deal  with  non-crop  activities  using  the  ME  methodol- 
ogy. As  a  result,  all  farms  involved  in  such  non-crop  activities  are  excluded  from  the  analy- 
sis. Along  the  same  line,  some  Saskatchewan  farms  were  single-enterprise  farms  and  were 
therefore  also  excluded. 
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A  common  approach  to  estimate  cost-allocation  coefficients5  from  whole-farm  accounting 
data  has  been  to  consider  a  system  of  linear  equations,  representing  the  derived  demand  for 
each  farm  input  as  a  function  of  several  farm  outputs.  In  other  words,  inputs  and  outputs, 
both  expressed  as  costs  and  revenues,  are  treated  as  the  dependent  and  independent  vari- 
ables, respectively.6 

Given  /  inputs  used  by  T  farms  to  produce  K  outputs,  the  system  of  linear  input-demand 
equations  can  be  written  as  follows  (see  also  Box  2  for  the  adopted  notation): 

K 

x't  =X%  ^  +u'<  fori  =  \X.,N  and  /  =  1,2,..T  (0 


k=\ 


where 


Xj  =      the  total  cost  associated  with  input  /  (i  =  1,  2,...,  I),  paid  by  farm  t, 

ailc  =  the  unknown  cost-allocation  coefficients,  which  are  defined  here  as  the 
average  (i.e.  for  all  farms)  expenditure  on  input  i  required  per  unit  of 
output  value  k, 

yk  =      the  total  revenue  corresponding  to  output  k  (k  =  1,  2,..,  K)  produced  by 

farm  t, 

Uj  =      a  random  disturbance,  which  is  specific  to  each  input  and  to  each  farm. 


5.  This  section  is  taken  in  part  from  Leon  et  al.  (1999),  where  the  use  of  the  ME  methodology  to  estimate  cost- 
allocation  coefficients  is  developed  at  full  length. 

6.  The  theoretical  foundations  of  such  a  linear  input  demand  model  rest  upon  the  adoption  of  a  Leontief 
(fixed  proportions)  technology  or  a  linear  production  process  framework  as  described  in  Ray  (1985)  and 
Scandizzo  (1990). 
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Box  2:  Notation 


„  Dimension 

Index  Domain  .       ..     ..  . 

(applications,  see  below) 

inputs 
i(l)  -excluding  fallow  costs  (application  1)  9 

-  including  fallow  costs  (application  2)  10 

k(K)  outputs  5 

t(T)  farms  30 

m(M)  support  values  for  the  ajk  parameters  1 1 

n(N)  support  values  for  the  errors  3 

Several  "classical"  tools  can  be  (and  were)  used  to  estimate  the  ajk  parameters  of  the  input- 
demand  system  in  (1),  such  as  linear  regression  techniques  (LRT),  bayesian  estimation 
techniques  (BES)  and  linear  programming  (LP).  However,  the  use  of  these  techniques  creates 
various  practical  problems.  For  example,  it  is  a  well-known  result  that  the  application  of  LRT 
may  lead  to  ajk  estimates  that  are  negative.  Also,  due  to  the  overall  accounting  constraint 
(which  equates  total  costs  to  total  revenues),  the  disturbance  terms  of  the  various  input- 
demand  equations  are  not  independent  from  each  other.  As  a  result,  the  system  of  input 
demand  equations  is  singular,  which  further  invalidates  the  use  of  the  LRT  technique 
(Bewley  1986).7 

To  ensure  the  non-negativity  of  the  estimated  production  coefficients,  one  may  consider 
using  constrained  estimation  procedures,  such  as  BES  or  inequality  restricted  least  squares 
methods  (see,  for  instance,  Moxey  and  Tiffin  1994).  These  alternative  methods,  however,  are 
very  cumbersome  to  implement  and  they  do  not  incorporate  the  overall  accounting  balance 
because  all  the  equations  in  the  system  are  treated  separately. 

The  approach  adopted  here  implies  the  reparameterization  of  the  linear  statistical  model  in 
(1)  for  the  ME  methodology,  so  that  all  the  non-negativity  and  adding-up  constraints  are 
incorporated  into  the  estimation  procedure.8 

To  do  so,  Section  3.1  gives  the  definition  of  the  entropy  measure  used  to  estimate  the 
statistical  input-demand  model  in  (1).  Then  the  standard  ME  and  ME-tobit  formulations  of 
the  cost-allocation  problem  are  presented  in  Sections  3.2  and  3.3,  respectively.  Finally, 
Section  3.4  presents  and  discusses  the  statistical  properties  of  the  ME  and  ME-tobit  systems 
of  equations. 


7.  System-of-equations  estimation  procedures,  such  as  the  iterative  seemingly-unrelated  estimation 
technique,  could  have  been  used  to  overcome  this  singularity  problem.  However,  such  an  estimation 
procedure  does  not  prevent  the  possibility  that  aik  estimates  can  be  negative. 

8.  The  ME  methodology  is  advocated  when  we  are  facing  "ill-posed"  problems  in  the  sense  of  insufficient 
sample  information  as  described  by  Golan  et  al.  (1996b).  More  specifically,  such  a  situation  arises  when  the 
number  of  parameters  exceeds  the  number  of  observations.  We  do  not  experience  this  problem  in  the  case 
of  our  sample  of  Saskatchewan  crop  farms. 
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3.1  Shannon's  entropy  measure 

Given  n  data  points,  the  entropy  measure,  introduced  by  Shannon  (1948),  is  expressed  in  terms  of 
unobserved  probabilities  p  =  [p1?  p2, ...,  pn]'  as 

H(p)=-XPJln(p,)  (2) 

where  p  Inp  =  0  for  p  =  0,  and  H(p)  reaches  a  maximum  when  p1  =  p2  =  ...  =  pn  =  1/n;  that  is, 
when  the  probabilities  are  uniform.  In  other  words,  maximizing  the  entropy  measure  H(p) 
amounts  to  choosing  the  probability  vector  p  that  is  closest  to  the  uniform  distribution  and 
yet  consistent  with  what  we  know,  that  is,  with  the  available  data  and  relevant  constraints 
(Jaynes  1957a,  b).9 

3.2  Standard  ME  formulation  of  the  cost-allocation  problem 

Expression  (1)  can  be  treated  either  as  a  single-equation  linear  statistical  model  or  as  a  sys- 
tem-of-equations  statistical  model. 

Single-equation  linear  statistical  model 

Given  the  above  entropy  measure  (2),  the  estimation  of  the  parameters  aik  and  the  error  process 
terms  ui  of  the  system  of  input-demand  equations  in  (1)  can  be  formulated  as  the  problem  of 
finding  the  expected  value  of  a  probability  distribution.  This  value  is  defined  over  a  set  of 
known  and  discrete  "support  values"  within  pre-specified  intervals  (Golan  et  al.  1996b). 
Thus,  for  each  input  /'  we  get  the  following  equations: 

M 

a*  =X-£  Pi  =-.  'A*,   for/  =  l,2,.../  wdk  =  \,2,..,K  (3) 

m=\ 


and 


N 


",'=XV>"=V"X>   fon  =  l,2,.../  and/  =  l,2,...,r  (4) 


«=i 


where 


Zjk         =      the  support  vectors  (of  dimension  M)  for  the  K  production  coefficients 
ajk  associated  with  input  /, 

vit         =      the  support  vectors  (of  dimension  /V)  for  the  Terror  terms  ui   associated 
with  input  i 

Pik.(wit)  =    tne  corresponding  unknown  probability  vectors  for  aitl  [«,-  ] . 


9.       It  can  easily  be  shown  that  the  implied  probability  distribution  is  the  one  that  "has  the  lowest  information 
content"  or  that  "can  be  realized  in  the  greatest  number  of  ways"  consistent  with  what  we  know. 
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The  elements  of  the  support  vectors  are  centred  about  the  likely  value  of  the  parameters  and 
noise  components  to  be  recovered.  For  our  purposes,  it  suffices  to  consider  common  sets  of 
support  values  for  the  unknown  production  coefficients  and  error  terms,  respectively;  that  is, 
Zik  =  z  =  [z1,  Z2,  ...,  ZM],  for  all  i,  k,  and  vit  =  v  =  [v1 ,  v2,  ...,  v^]  for  all  i,  t.  Doing  so  simplifies  our 
notation  considerably,  so  that  the  expressions  (3)  and  (4)  can  be  rewritten  as  follows: 


M 


*,*=!> >;=-->,*,    for  i  =  1,2,.... K 


(5) 


and 


N 


u\  =£v"  w{  =v'w„,  for  i  =  1,  2, /and  k  =1,2, 


K 


(6) 


n=\ 


System-of-eq nations  statistical  model 

Alternatively,  expression  (1)  can  be  treated  as  a  system  of  interdependent  (seemingly-unre- 
lated) equations.  Doing  so  means  that  all  the  inputs  are  taken  into  account  simultaneously10 
Hence,  an  additional  accounting  restriction  should  be  imposed  for  each  type  of  output  k: 


^alk=\,    for  all  k 


(7) 


i=i 


This  restriction  ensures   that  the  overall  adding-up  consistency  or  accounting  balance 
between  total  revenue  and  total  cost  is  always  satisfied-as  explained  below. 

Given  the  accounting  balance  between  total  revenue  and  total  cost  for  each  farm,  a  singular- 
ity problem  occurs.  Specifically,  given  the  /  equations  in  (1)  for  each  farm,  we  can  write: 


/  /        (   K  \ 


1=1 


A-=l 


(8) 


/ 


Introducing  the  cross-equation  restriction  E.-a.-i,  =   1  implies  that 

i*:'=i  ft  «.<  V;+x  «;=i>;+i «.' 

i=l  k=\     y  i=l  J  i=l  k=\  1=1 

Since,  by  construction, 


(9) 


X*;=I  yl 


(10) 


1=1 


A=l 


10.    Considering  the  /  linear  models  in  (1)  simultaneously,  for  each  farm  t,  implies  that  we  assume  that  the  errors  u\  are 
"contemporaneously"  correlated  (i.e.  for  each  individual  farm),  but  uncorrelated  across  the  farms. 
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this  implies  that 


(ii) 


;=1 


Hence,  the  system  is  singular  and  the  inverse  of  the  variance-covariance  matrix  does  not 
exist.  Consequently,  ordinary  least  squares  (OLS)  or  maximum  likelihood  (ML)  procedures 
cannot  be  used  properly  to  estimate  the  unknown  parameters  of  the  system  of  equations. 

With  the  above  specifications,  the  following  ME  problem  can  now  be  formulated: 

given  the  support  vectors  [z,  v]  along  with  the  sample  data  on  inputs  xt  and  outputs  yk  , 

find  the  probability  vectors/?  »  0  and  w  »  0  that  maximize  the  following  entropy  measure: 


H(p,  w)  =  -p1  In  p  -  w'  In  w 


(12-1) 


subject  to  the  data-consistency  constraints, 
x  =  Ya  +  u  =  Y(z'p)  +  v'w 


(12-2) 


the  adding-up  constraints, 

M 

£rf=l,   for  all  uk 


m=\ 


(12-3) 


N 


<=1,  for  all/,/ 


(12-4) 


and  the  additional  cross-equation  or  accounting  restriction 
IX  =  Xr,#*=1>  for  all /c 


7=1 


(=1 


(12-5) 


The  solution  to  the  above  ME  optimization  problem  is  unique  and  the  unknown  probability 
vectors  p  and  w  can  be  derived  by  using  non-linear  programming  software,  such  as  the 
computer  program  GAMS  (General  Algebraic  Modelling  System).11 

3.3     ME-tobit  formulation  of  the  cost-allocation  problem 

If  some  cost  items  are  equal  to  zero,  as  could  be  observed  for  a  number  of  farms,  the  use  of  a 
censored-data  or  tobit  variant  of  the  linear  statistical  model  is  recommended,  where  the 
observations  are  ordered  as  follows: 


x  =  Y  A  +  u  = 


x,>0 
x2=0 

= 

Y2_ 

A  + 

_U2. 

(13) 


11.    For  more  details  on  this  computer  program,  see  Brooke  et  al.  (1988). 
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where  x{  is  observed  and  x2  is  unobserved,  and  Y^  and  Y2  are  the  corresponding  output 
matrices.  Using  the  ME  formalism  yields  the  following  reparameterized  model  (Golan  et  al. 
1996b,  p.  270): 


x,>0 
x2=0 


LY2 


z'p  + 


V,  w, 


v2  w2 


(14) 


which  leads  to  the  following  ME-tobit  problem: 

given  the  support  vectors  [z,  wv  v2]  along  with  the  sample  data  on  inputs  xi  and  outputs 
yk  ,  find  the  probability  vectors  p  »  0,  Wj  »  0,  and  w2  »  0  that  maximize  the  following 
entropy  measure: 


H(p,w,,w2)  =  -p,lnp-w1'lnwl  -w2'lnw. 


(15-1) 


subject  to  the  data-consistency  constraints, 
x,  =Y1z'p  +  v1'w1 


(15-2) 


0>Y2z*p  +  v2,w 


2     vv  2 


(15-3) 


and  the  adding-up  constraints, 


M 

2_,Pik=^     for  all  i,  A: 


m=l 


(15-4) 


I 


w".t=l,    for  all  /and/  =1,  2,  ...J 


(15-5) 


I 
2_.w2it  =  lj   for a11  *>  and  f  =  7'  2>—T2>  where  T]  +  T2=T 


n=\ 


(15-6) 


and  the  accounting  constraint 


I*  ■  I-- 


A*=l 


/=i 


(15-7) 


Given  the  above  specification,  there  are  T{  (farms  with)  positive  observations  and  T2  (farms 
with)  zero  observations  for  input  i. 
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Again,  the  above  ME  optimization  problem  can  be  solved  along  similar  lines  as  the  problems 
previously  defined,  using  standard  non-linear  programming  software,  such  as  GAMS. 

3.4     Statistical  properties  of  the  ME  and  ME-tobit  systems  of 
equations 

The  statistical  properties  of  ME  and  ME-tobit  systems  of  equations  are  assessed  by  three  indi- 
cators: the  entropy  ratio  statistic,  the  normalized  entropy  indicator  and  asymptotic  standard 
errors. 

Entropy  ratio  statistic 

The  statistical  properties  of  the  ME  estimators  of  the  (constrained  or  unconstrained)  system 
of  equations  were  recently  investigated  by  Golan  and  Judge  (1996a),  Golan  et  al.  (1996b, 
2001),  Marsh  et  al.  (1998)  and  Mittelhammer  and  Cardell  (1997).  Assuming  mild  conditions 
on  the  error  terms,  on  the  support  values  z  and  v,  and  on  the  matrix  of  explanatory  variables 
Y,12  it  can  be  shown  that  the  ME  estimators  of  the  aik  parameters  are  consistent  and  asymp- 
totically normal.  Given  this  result,  it  is  possible  to  show  that  the  entropy  ratio  statistic  for  the 
different  parameters  of  the  unknown  distribution  generating  the  data  distributions  has  a  lim- 
iting distribution.  Applied  to  the  system  of  linear  input-demand  equations,  the  entropy  ratio 
(ER)  statistic  is  defined  as  follows. 

Under  the  null-hypothesis  "//0  =  the  sum  of  the  output  coefficients  ajk  ,  for  a  given  output  k, 
is  equal  to  one  {^kaik  =  1 )  ",  the  entropy  ratio13  is  defined  as: 

ER{lka^=1}=2SuN(lkalk^l)-2SR(Xka^=l)  (16) 

where  SUN  and  SR  are  the  values  of  the  unrestricted  and  restricted  ME  functions,  respectively, 
as  given  by  the  expressions  (12-1)  and  (15-1),  respectively.  Under  the  null-hypothesis  //0,  the 
entropy  ratio  (ER)  statistic  follows  a  x2  distribution,  with  K  degrees  of  freedom,  where  K  is 
the  number  of  constraints  imposed  on  the  aik  parameters. 

Normalized  entropy  indicators 

The  relative  performance  of  the  ME  estimation  of  the  ajk  coefficients  can  be  measured  by  the 
normalized  entropy  indicator,  defined  as  the  proportion  of  the  total  remaining  uncertainty. 
The  normalized  entropy  measure,  for  I*K  parameters  (/  inputs  and  K  outputs)  and  M  possi- 
ble outcomes  (M  support  values  for  the  ajk  parameters),  is  defined  by  Golan  et  al.  (1996b, 
p.  93)  as  follows: 

S(p)=^M-  (17) 

IKln(M) 


12.  Specifically,  there  are  four  conditions  (Golan  et  al.  2001,  p.  24):  the  error  support  values  v  are  symmetric 
around  zero,  the  support  values  z  associated  with  the  aik  parameters  span  the  true  value  for  each  of  the 
unknown  parameters  and  have  finite  lower  and  upper  bounds,  the  errors  are  independently  and  identically 
distributed  with  mean  zero  and  with  a  contemporaneous  N  x  N  variance-covariance  matrix  I ;  and  the 
matrix  Y'Y  exists  and  is  non-singular. 

13.  It  is  worth  noting  the  similarity  between  this  entropy  ratio  statistic  and  the  likelihood  ratio  statistic. 
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where  £(/?)£ [0,  1],  and  where  IKln(M)  represents  the  maximum  uncertainty  —  that  is,  the 
entropy  level  of  the  uniform  distribution.  When  S(p)  =  0,  there  is  no  uncertainty;  when 
S(p)  =  1 ,  the  distribution  of  the  K>7  parameters  is  identical  to  the  uniform  prior  distribu- 
tion. Along  similar  lines,  the  informational  content  of  the  noise  or  error  component  can  be 
assessed  through  the  normalized  entropy  measure,  for  /*T  errors  (J  inputs  and  T  farms)  and 
N  possible  outcomes  (N  support  values  for  the  error  terms),  defined  by  Golan  et  al.  (1996b,  p. 
93)  as  follows: 

S(w)=^^  (18) 

ITln(N) 

where  5(w)e[0,  1],  and  where  ITln(N)  represents  the  maximum  uncertainty-that  is,  the 
entropy  level  of  the  uniform  distribution. 

Asymptotic  standard  errors 

The  asymptotic  standard  errors  of  the  estimated  parameters  ajk  are  calculated  by  defining, 
in  the  first  stage,  the  variance-covariance  matrix  of  the  estimated  residuals,  denoted  by  E . 
The  elements  of  £  are  defined  as  follows:14 


o,=-2,u,tuij. 
a  t 


(19) 

.       1T 

where 


a*-Xw2v;  (20) 

n 

Given  the  estimated  variance-covariance  matrix   Z,  the  asymptotic  variance-covariance 
matrix  of  the  constrained  parameters  aik  is: 


^r-gme  =«-nR,(Ri2R)",Ri2  (21) 


where 


Q  =  plimT-1[Y'(£-1  ®I)y]"'  (22) 

The  matrix  /?,  of  dimension  K*I,  is  defined  in  such  a  way  (containing  "l's"  in  the  appropriate 
positions)  that  Ra  =  e,  where  a  is  the  vector  of  aik  parameters  and  e  is  a  unit  vector.  As  a 
result,  the  complete  set  of  K  equality  conditions  that  the  sum  of  the  aik  parameters  equal  one, 
for  each  output  k  (k=  1,.. .,K),  are  satisfied. 


14.     It  is  a  well-known  result  that  for  small  samples  the  use  of  1/T  in  expression  (19)  yields  a  biased  estimate  of 
cr  .  To  remedy  this  problem,  Zellner  and  Huang  (1962)  suggested  the  replacement  of  T  by  (T-Kf  (T-Kj) 
in  expression  (19),  where  AT,  and  Kj  are  the  number  of  explanatory  variables. 
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Next,  it  can  be  established  that  the  restricted  ME  estimator  of  the  aik  parameters  is  consistent 
and  asymptotically  normally  distributed  (Golan  et  al.  2001): 

VT(a-a)-^N(0,QR_OME)  (23) 

All  these  statistical  results  are  used  in  applying  the  ME  methodology  to  Saskatchewan  farms 
and  are  reported  in  Section  4. 
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Section  4:   Maximum  entropy  applied  to  the 

Saskatchewan  data  sample 

The  ME  methodology  is  applied  to  the  accounting  data  for  30  crop  farms  in  Saskatchewan.6 
As  some  farms  in  the  sample  exhibit  zero  and/ or  negative  values  for  certain  input 
categories/  the  standard  ME  methodology  had  to  be  adjusted.  Specifically,  the  tobit 
formulation  of  the  ME  methodology  (ME-tobit),  as  given  by  expression  (13),  is  used  and 
applied  to  the  data  sample  of  farms.8 

The  data  are  for  nine  input  categories  (7  =  9)  and  five  output  categories  (K  =  5).  All  input 
expenditures  (accounting  costs)  and  output  values  are  expressed  in  Canadian  dollars  ($). 

This  section  has  six  sub-sections.  Section  4.1  has  three  alternative  ME-tobit  model  specifica- 
tions with  varying  support  sets.  In  Section  4.2,  these  alternative  model  specifications  are  then 
applied  to  two  variants  of  the  Saskatchewan  data  sample  — one  excluding  fallow  costs  as  a 
separate  cost  item,  the  other  including  the  fallow  costs  as  a  separate  cost  item.  We  present  the 
estimation  results  of  the  various  ME-tobit  model  specifications  in  Section  4.3.  We  discuss  and 
validate  the  performance  of  these  alternative  ME-tobit  model  specifications  in  the  last  three 
sub-sections. 

4.1     Three  alternative  model  specifications  with  varying  support 
sets 

Three  alternative  model  specifications  are  used  with  varying  designs  of  the  support  sets  for 
the  cost-allocation  coefficients.  The  three  sets  of  support  values  (intervals)  are  different  only 
with  regard  to  the  width  of  the  interval  and/ or  the  type  of  spacing  of  the  support  values 
within  the  chosen  interval.  An  overview  of  the  three  model  specifications,  designated  here 
as  Models  A,  B  and  C  and  their  associated  support  sets  are  given  in  Table  2.  For  each  model, 
eleven  values  (M  =  11)  were  chosen  for  the  common  parameter  support  vector  z  (for  all  i,  k), 

6.  Given  the  small  size  of  the  sample  (only  30  observations),  the  cr  's  are  corrected  for  the  number  of  degrees 
of  freedom  using  Zellner  and  Huang's  suggestion  (see  Footnote  14). 

7.  This  situation  occurs  for  the  four  input  categories:  fertilizers  (one  farm),  fuel,  power  and  grease  (one  farm), 
paid  salaries  (six  farms)  and  net  operating  income  (three  farms). 

8.  The  standard  ME  applications  were  also  conducted  but  the  results  which  are  not  reported  here  are  available 
from  the  authors  upon  request.  For  other  recent  applications  of  the  ME-tobit  approach  to  systems  of 
equations  in  the  context  of  agricultural  economics,  see  Arndt  (1999)  and  Golan  et  al.  (2001) 
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along  with  three  values  (N  =  3)  for  the  common  error  support  vector  v  (for  all  i,  t).  The 
selected  values  of  the  common  support  vector  z  are  motivated  by  "prior  information"  about 
the  range  of  magnitudes  for  the  cost-allocation  coefficients.10 

In  Model  A,  the  interval  is  defined  in  such  a  way  that  all  the  support  values  are  within  the 
admissible  [0.0,  1.0]  range,  while  being  equally  spaced  over  the  interval.  Model  B  is 
grounded  on  the  plausible  assumption  that  the  coefficients  should  be  closer  to  zero  than  to 
one.  Hence,  the  size  of  the  interval  is  reduced,  spanning  the  range  [0,  0.5].  In  addition,  the 
support  values  are  equally  spaced  over  the  reduced  interval.  In  Model  C,  a  so-called  "asym- 
metric" and  "left-skewed"  type  of  spacing  over  the  reduced  interval  [0,  0.5]  is  considered, 
with  a  larger  number  of  support  values  that  are  close  to  zero  and  fewer  support  values  that 
are  close  to  0.5. 

The  error  support  vector  v  is  treated  differently,  depending  on  whether  all  the  observed  val- 
ues of  the  input  demand  are  positive  or  whether  some  of  them  are  negative  or  zero.  In  the 
first  case  (all  positive),  the  selected  support  values  are  based  on  the  results  of  prior  OLS  esti- 
mations. The  support  values  are  always  symmetrically  defined  around  zero.  The  endpoints 
[-  3OM0,  +  3ou]  of  the  intervals  are  based  on  the  so-called  "three-sigma  rule",  where  the 
range  was  derived  from  prior  OLS  standard-error-of-regression  results.11  In  the  second  case 
(some  negative  or  zero),  the  OLS  estimate  of  the  standard-error-of-regression  is  biased  and 
hence,  cannot  be  used  to  generate  the  support  range  in  the  ME-tobit  model.  Therefore,  an 
alternative  estimate  proposed  by  Golan  et  al.  (1997)  is  used,  which  is  based  on  an  assumed 
uniform  distribution  of  the  data.12, 13 

4.2     Two  applications  excluding  and  including  fallow  costs 

Application  1:  TJie  ME-tobit  model  excluding  fallow  costs  as  a  separate  cost  item 

The  first  application  implies  a  tobit  version  of  the  ME  methodology  (ME-tobit),  where  input 
expenditures  are  linked  to  five  crop  outputs,  expressed  in  monetary  terms.  This  application 
leads  to  the  Models  Al,  Bl  and  CI.  Fallow  costs  as  such  are  included  in  the  estimation  (but 
not  as  a  separate  cost  item).  In  all  three  model  variants,  it  is  implicitly  assumed  that  all  the 


9.  For  a  more  thorough  analysis  on  the  sensitivity  of  the  ME  estimates  to  the  width  of  the  support  interval 
(end  points)  and/ or  the  type  of  spacing  of  the  support  values  within  the  pre-specified  intervals,  see  Leon  et 
al.  (1999,  p.  435-437). 

10.  The  "prior  information"  relates  to  the  fact  that  all  the  cost-allocation  coefficients  must  lie  within  the  range 
[0, 1].  Furthermore,  common  sense  and  the  existence  of  a  large  number  of  cost  items  (nine  or  ten  depending 
upon  the  applications)  lead  us  to  conclude  that  the  estimated  aik's  should  be  closer  to  zero  than  one.  This 
fact  is  confirmed  by  OLS  estimation  experiments  that  all  estimated  coefficients  were  close  to  zero  and 
smaller  than  0.5.  All  this  factual  evidence  also  suggests  that  using  left-skewed  distribution  for  the  values  of 
the  common  support  vector  z  is  well  justified. 

11.  Concerning  the  use  of  this  "three-sigma  rule,"  we  followed  the  recommendations  of  Golan  et  al.  (1996b,  p. 
88)  who  used  it  in  their  econometric  experiments. 

12.  This  alternative  estimate  of  cu  is  obtained  by  assuming  a  uniform  and  censored  distribution  of  the  data. 
Let  xmax  denote  the  largest  value  of  xt .  Then,  knowing  the  range  in  the  sample  and  the  proportions  of  the 
censored  observations  (that  is,  where  x~  0)  in  the  data  sample,  the  unobserved  xmin  can  be  derived.  The 
uniform  variance  sx  =  (xmax  -  xmin )    /12  is  used  as  an  estimator  of  the  variance  au2  . 

13.  In  terms  of  software  used,  the  system  of  input-demand  equations  was  estimated  using  the  GAMS  computer 
program,  while  the  variance-covariance  matrix  of  the  restricted  ME  parameters  was  obtained  using  the 
econometric  Time  Sharing  Processor  (TSP)  computer  package  (Hall  and  Cummins  1998).  Examples  of  these 
two  computer  programs  are  provided  in  Appendices  C  and  D,  respectively. 
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Model 
specification 

Number  of 
support  values 

Type  of 
spacing 

Model  A 
Model  B 

11 
11 

Symmetric 
Symmetric 

inputs  (costs)  associated  with  fallowed  land  operations  are  imputed  to  each  crop  output,  and 
are  invariant  with  respect  to  time  — that  is,  there  is  no  year-to-year  variation.  Given  this 
assumption,  there  is  no  need  to  include  the  number  of  hectares  of  fallowed  land  as  an 
additional  explanatory  variable  in  the  problem  formulation.  Ignoring  fallow  costs  as  a 
separate  cost  item  is  not  entirely  satisfactory  because  the  resulting  cost-allocation  coefficients 
are  considerably  overestimating  their  true  values. 

Table  2:    Model  specifications  and  corresponding  support  sets  for  the 
unknown  cost-allocation  parameters  aik 

Selected  support  values 

0.0.1,0.2,0.3,0.4.0.5,0.6,0.7,0.8,0.9,1  0 

0,0  05,0,10,0.15,0.20,0  25,0.30.0  35,0  40,0.45,0  50 

Non-symmetric  left-         0,0.025,0  050,0.075,  0.100,  0.125,  0,150,  0.200,  0  300, 
Model  C  11  skewed  0.400,0  500 

Application  2  :  The  ME-tobit  model  including  fallow  costs  as  a  separate  cost  item 

The  second  application  implies  a  tobit  version  of  the  ME  methodology  (ME-tobit)  where,  in 
contrast  with  the  first  application,  fallow  costs  are  included  as  a  separate  cost  item.  Specifi- 
cally, an  additional  input  category,  called  "fallow  costs,"  is  incorporated  into  the  model.14 
This  input  category  depends  on  the  five  categories  of  crop  outputs.15  As  a  result,  the  model 
contains  ten  input-demand  equations  and  the  three  estimated  model  specifications  are  now 
designated  as  Models  A2,  B2  and  C2.  For  purposes  of  comparison,  the  ME-tobit  model  was 
also  estimated  for  the  three  different  sets  of  support  values  defined  earlier. 

4.3     Results  of  the  ME-tobit  estimation 

The  estimation  results  are  presented  in  Tables  3  and  4  for  the  ME-tobit  Applications  1  and  2, 
respectively.  Additional  results  on  the  statistical  tests  and  the  performance  of  the  models  are 
reported  in  Tables  5,  6,  7  and  8.  The  empirical  results  are  quite  reasonable  in  all  respects, 
regardless  of  the  model  specification  used.  The  only  issue  of  concern  is  the  relatively  weak 
statistical  significance  of  the  estimated  parameters  (see  Table  5).  Looking  at  Model  Al,  it  can 
be  seen  that  26  out  of  45  estimated  coefficients  (nine  inputs  times  five  outputs)  are  statisti- 
cally significant. 


14.  Six  farms  out  of  30  do  not  use  summerfallow  and  hence  are  characterized  by  fallow  costs  equal  to  zero. 

15.  The  costs  associated  with  fallow  operations  were  obtained  using  the  cropping  enterprise  database  which 
records  all  the  expenses  imputed  to  the  maintenance  of  fallowed  land.  These  operations  concern  mainly 
pesticides,  some  other  direct  material  inputs,  fuel,  power  and  grease  and  repairs.  In  the  category  of  paid 
salaries,  the  imputation  of  the  cost  attributable  to  fallowed  land  was  undertaken  by  taking  the  proportion  of 
management  fees  devoted  to  fallow  land  appearing  in  the  cropping  enterprise  database.  Finally,  in  other 
fixed  cash  expenses,  cash  cost  associated  with  fallowed  land  was  measured  as  follows:  interest  paid, 
property  taxes,  overhead  materials,  and  cash  land  lease  that  can  be  imputed  to  fallow  maintenance 
operations  were  computed  for  each  farm  using  the  cropping  enterprise  database.  These  expenses  were  then 
converted  into  percentages  which  were  applied  to  total  fixed  cash  expenses.  Finally,  to  ensure  that  the 
accounting  identity  between  total  revenues  and  costs  is  satisfied  for  each  farm,  all  expenses  other  than 
fallow  costs  were  adjusted  downward  by  an  amount  imputed  to  fallow  costs. 
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Specifically,  14  coefficients  are  significant  at  the  5  percent  level  of  significance,  while  12  coef- 
ficients are  significant  at  the  10  percent  and  15  percent  levels  of  significance,  respectively. 
Examining  each  input-demand  equation  separately,  it  can  be  seen  that  three  to  four  out  of 
five  coefficients  are  significant  in  the  equations  for  seeds,  pesticides,  fuel,  power  and  grease 
and  other  fixed  cash  expenses,  whereas  only  one  coefficient  is  significant  in  the  equations  for 
fertilizers,  other  direct  material  inputs  and  repairs.  In  the  equation  for  net  operating  income, 
only  one  coefficient  is  significant  at  the  10  percent  level  of  significance.  A  similar  pattern  is 
discernible  for  the  other  model  specifications. 

Consistent  estimation  of  the  coefficients  was  possible  by  imposing  the  adding-up  constraint, 
which  also  ensures  that  the  sum  of  the  coefficients,  for  each  output  separately,  is  equal  to 
one.  These  restrictions  are  tested  by  comparing  the  restricted  and  the  unrestricted  models, 
using  the  entropy  ratio  (ER)  statistic.  The  null  hypothesis,  stating  that  the  sum  over  all  inputs 
of  the  ajk  coefficients  is  equal  to  one,  cannot  be  rejected  if  the  calculated  statistic  ER  is 
smaller  than  11.09  — which  is  the  critical  x  value  at  the  5  percent  level  of  significance.  An 
inspection  of  these  test  results,  provided  in  Table  6,  reveals  that  the  null-hypothesis  can  be 
rejected  only  for  Models  Al  and  A2. 

The  normalized  entropy  (NE)  indicators  are  important  criteria  for  assessing  the  informa- 
tional content  of  each  model  specification.  Specifically,  a  model  is  viewed  as  "superior"  if  its 
associated  NE  measure  S(p)  is  greater  than  the  values  of  the  same  indicator  obtained  for 
alternative  model  specifications.  This,  means  that  the  "superior"  model  would  yield  a  solu- 
tion for  the  recovered  cost-allocation  coefficients  that  is  more  "uniform"  or  closer  to  the  prior 
distribution.  Using  this  criterion,  S(p)  increases  from  0.55  in  Model  Al  to  0.77  in  Model  CI 
(see  Table  7).  Hence,  we  can  conclude  that  a  "left-skewed"  spacing  of  the  support  values  over 
a  reduced  interval  [0,  0.5]  yields  more  satisfactory  results.  This  finding  is  not  surprising  how- 
ever, and  is  in  line  with  our  prior  judgement  that  the  coefficients  are  more  likely  to  be  closer 
to  zero  than  to  one.  A  computation  of  the  normalized  entropy  indicator16  for  all  the  recov- 
ered cost  allocation  coefficients  (the  ciik's)  confirms  this  point  (see  Appendix  E  for  more 
details  on  results).  A  similar  pattern  holds  in  the  case  of  Application  2  where  Model  C2  is 
preferred  to  Model  A2. 

Comparing  the  two  "best  performing"  model  specifications,  that  is,  Models  CI  and  C2,  then 
Model  C2  is  "superior"  because  it  shows  a  higher  value  (0.7706  against  0.8363)  (see  Table  7) 
for  the  NE  measure  S{p) .  With  regard  to  the  informational  content  of  the  noise  ratio,  the  NE 
measure  S{w)  is  invariant  to  the  choice  of  the  support  set  (see  Table  7).  This  result  is  not  sur- 
prising, as  the  initial  support  values  for  w  have  not  been  changed  across  the  various  model 
specifications. 


16.    The  normalized  entropy  ratio  associated  to  each  ailc  is  defined  as  follows:  L,r* 

SAP)=     m\n(M) 
where  M  is  the  number  of  support  values  equal  to  eleven. 

The  interpretation  of  Sik(p)  is  the  same  as  for  S(p) :  if  it  is  equal  to  one,  it  means  that  the  distribution  of  the 
coefficient  aik  is  identical  to  the  uniform  prior  distribution.  An  examination  of  the  Sik(j))  's  in  Appendix  E 
shows  that  all  of  them  increase  significantly  when  we  switch  from  Models  A  to  B.  This  finding  confirms  the 
fact  that  shrinking  the  range  of  the  support  values  leads  to  cost  allocation  coefficients  which  are  more  in 
line  with  the  prior  belief  that  the  aik  's  are  closer  to  zero  than  one.  This  pattern  is  also  reinforced,  but  to  a 
lesser  extent,  when  we  move  from  Model  specifications  B  to  C.  However,  this  is  not  valid  for  all  the 
estimated  coefficients  aik  .  In  fact,  few  of  them  have  their  normalized  entropy  indicator  declining  in  value. 
This  latter  pattern  indicates  that  the  symmetric  type  of  spacing  of  support  values  is  more  appropriate  for 
those  cost  allocation  coefficients. 
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Table  3:   Application  1 — ME-tobit  estimates  of  cost-allocation  coefficients 
(excluding  fallow  costs  as  a  separate  cost  item) 


Model  A1 

Model  B1 

Model  C1 

Coefficients 

Standard 
errors 

Coefficients 

Standard 
errors 

Coefficients 

Standard 
errors 

1.  Seeds 

Wheat 

0067 

00218 

0062 

00219 

0063 

0.0222 

Other  grains 

0013 

00307 

0020 

0  0297 

0020 

0.0313 

Canola 

0034 

00183 

0034 

00177 

0034 

00187 

Other  oilseeds 

0  069 

0.0645 

0.075 

0.0629 

0072 

0.0663 

Other  crops 

0  066 

00157 

0069 

00152 

0  067 

00160 

2.  Fertilizers 

Wheat 

0  194 

0  1237 

0.192 

00481 

0  193 

0.1276 

Other  grains 

0.158 

01791 

0153 

0.0676 

0  156 

01795 

Canola 

0  204 

0  1072 

0.206 

0.0404 

0.206 

0  1074 

Other  oilseeds 

0  090 

0.3796 

0094 

0.1433 

0.090 

0  3806 

Other  crops 

0.112 

0.0917 

0.117 

0  0346 

0.114 

0.0920 

3.  Pesticides 

Wheat 

0088 

0.0430 

'  0.101 

0.0381 

0.090 

0.0439 

Other  grains 

0005 

0.0604 

0006 

0.0535 

0010 

0.0617 

Canola 

0.057 

0  0362 

0.053 

0  0320 

0.056 

0.0369 

Other  oilseeds 

0400 

01281 

0.356 

0.1135 

0.386 

01309 

Other  crops 

0096 

0.0307 

0091 

0  0274 

0.095 

0.0316 

4.  Other  direct  material  inputs 

Wheat 

0034 

0.0215 

0032 

00189 

0.033 

0.0217 

Other  grains 

0.047 

0.0302 

0.051 

0  0266 

0.050 

0.0306 

Canola 

0.029 

0.0180 

0.029 

0.0159 

0028 

0.0183 

Other  oilseeds 

0  032 

0.0641 

0.041 

0.0563 

0040 

0.0648 

Other  crops 

0061 

0.0155 

0063 

00136 

0  062 

00157 

5.  Fuel,  power  and  grease 

Wheat 

0.055 

00305 

0057 

00315 

0.059 

0.0343 

Other  grains 

0114 

0  0428 

0  124 

0  0443 

0117 

0.0482 

Canola 

0062 

00256 

0.070 

0  0265 

0  065 

00289 

Other  oilseeds 

0071 

0  0908 

0.077 

0.0939 

0.073 

0  1023 

Other  crops 

0.063 

0.0219 

0.075 

0.0227 

0.070 

0.0247 

6.  Repairs 

Wheat 

0.089 

00312 

0.080 

0.0285 

0081 

00321 

Other  grains 

0071 

0  0438 

0.080 

0  0400 

0.076 

0.0451 

Canola 

0017 

0.0262 

0023 

0.0240 

0.023 

00270 

Other  oilseeds 

0076 

0.0929 

0.080 

0.0849 

0076 

0.0957 

Other  crops 

0  034 

00225 

0.040 

0.0205 

0  039 

0  0231 

7.  Paid  salaries 

Wheat 

0077 

0  1080 

0  091 

00909 

0  088 

01149 

Other  grains 

0.170 

01519 

0169 

0  1278 

0166 

01617 

Canola 

0116 

0.0909 

0143 

0.0765 

0125 

0.0967 

Other  oilseeds 

0  064 

0.3220 

0074 

02709 

0.071 

0  3427 

Other  crops 

0  135 

00778 

0157 

0.0655 

0143 

00828 

8.  Other  fixed  cash  expenses 

Wheat 

0  196 

0.0610 

0.180 

0  0520 

0  190 

0.0618 

Other  grains 

0211 

00857 

0.199 

00731 

0204 

00870 

Canola 

0125 

0.0513 

0141 

0.0438 

0131 

0  0521 

Other  oilseeds 

0  120 

01818 

0115 

01551 

0110 

0  1844 

Other  crops 

0.101 

0  0439 

0114 

00375 

0  107 

0  0446 

9.  Net  operating  income 

Wheat 

0  199 

02503 

0.204 

0  1854 

0202 

0  2660 

Other  grains 

0211 

0  3520 

0198 

0.2608 

0.200 

03741 

Canola 

0356 

02106 

0301 

0.1561 

0  332 

0.2239 

Other  oilseeds 

0077 

0  7462 

0  088 

05529 

0081 

07931 

Other  crops 

0  332 

0  1803 

0273 

0  1336 

0  303 

0  1917 
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Table  4:      Application  2 — ME-tobit  estimates  of  cost-allocation  coefficients 
{including  fallow  costs  as  a  separate  cost  item) 


Model  A2 

Model  B2 

Model  C2 

Coefficients 

Standard 
errors 

Coefficients 

Standard 
errors 

Coefficients 

Standard 
errors 

1 . Seeds 

Wheat 

0065 

00222 

0.061 

0  0226 

0.058 

0.0230 

Other  grains 

0015 

00313 

0022 

00318 

0.028 

0.0324 

Canola 

0.035 

00187 

0035 

00190 

0.035 

00194 

Other  oilseeds 

0  069 

0.0663 

0.073 

0  0674 

0071 

0.0688 

Other  crops 

0063 

00160 

0.066 

0.0163 

0.065 

0.0166 

2.  Fertilizers 

Wheat 

0  086 

0.0492 

0097 

00599 

0.099 

0.0579 

Other  grains 

0.125 

0.0693 

0  130 

00843 

0  129 

0.0815 

Canola 

0.1270 

.0414 

0141 

0.0504 

0135 

0  0487 

Other  oilseeds 

0  069 

0  1468 

0.075 

0  1787 

0073 

0.1728 

Other  crops 

0  102 

0.0355 

0115 

0.0432 

0.112 

00418 

3.  Pesticides 

Wheat 

0069 

00425 

0.070 

0  0436 

0  064 

0  0434 

Other  grains 

0002 

00597 

0.002 

00614 

0.006 

0.0616 

Canola 

0.063 

00357 

0065 

0  0367 

0  065 

0.0368 

Other  oilseeds 

0  385 

0  1266 

0.352 

0  1301 

0  365 

0.1305 

Other  crops 

0.104 

00306 

0  105 

0.0314 

0.108 

0.0315 

4.  Other  direct  material  inputs 

Wheat 

0036 

00213 

0.035 

00216 

0035 

0.0218 

Other  grains 

0045 

00300 

0.049 

0  0304 

0050 

0.0307 

Canola 

0  028 

00180 

0029 

00182 

0029 

0.0184 

Other  oilseeds 

0  030 

0.0636 

0.038 

0  0645 

0.044 

0.0650 

Other  crops 

0061 

0  0154 

0.062 

00156 

0061 

0  1572 

5.  Fuel,  power  and  grease 

Wheat 

0  060 

00301 

0064 

0  0384 

0067 

00389 

Other  grains 

0  106 

00424 

0.112 

0  0540 

0.110 

0.0453 

Canola 

0  060 

00254 

0.070 

0.0323 

0  069 

0.0271 

Other  oilseeds 

0  068 

0.0899 

0.071 

0.1144 

0.070 

0.0960 

Other  crops 

0  059 

00217 

0070 

00277 

0071 

0.0232 

6.  Repairs 

Wheat 

0  089 

0.0304 

0  082 

00314 

0.076 

00322 

Other  grains 

0.064 

0  0428 

0  072 

0  0442 

0073 

00453 

Canola 

0015 

0.0256 

0021 

0  0265 

0.026 

00271 

Other  oilseeds 

0  076 

0.0907 

0077 

00938 

0073 

0.0960 

Other  crops 

0032 

00219 

0037 

0.0227 

0041 

0.0232 

7.  Paid  salaries 

Wheat 

0.090 

0  1085 

0.110 

0.1282 

0.110 

0.1241 

Other  grains 

0.162 

00428 

0  157 

0  1804 

0.155 

0.1746 

Canola 

0  119 

00256 

0.150 

0  1079 

0.142 

0.1045 

Other  oilseeds 

0061 

00907 

0069 

0  3823 

0  069 

0.3701 

Other  crops 

0  127 

00219 

0143 

0  0924 

0  138 

0.0895 

8.  Other  fixed  cash  expenses 

Wheat 

0  193 

01085 

0182 

0  0520 

0  185 

0.0520 

Other  grains 

0.176 

0.1526 

0.169 

0  0732 

0.168 

00731 

Canola 

0.116 

00913 

0.130 

00438 

0  128 

00437 

Other  oilseeds 

0117 

0.3235 

0  108 

0.1551 

0.099 

0.1550 

Other  crops 

0  094 

00782 

0  103 

0  0375 

0  102 

00375 

9.  Fallow  cost 

Wheat 

0023 

0.0507 

0.025 

0  0274 

0  028 

0.0294 

Other  grains 

0  096 

00713 

0.098 

0  0386 

0092 

00414 

Canola 

0018 

00426 

0.022 

0  0231 

0025 

0  0248 

Other  oilseeds 

0  049 

01511 

0053 

00817 

0054 

00877 

Other  crops 

0024 

00365 

0030 

00197 

0033 

00212 

10.  Net  operating  income 

Wheat 

0  300 

02069 

0273 

0  2532 

0278 

02494 

Other  grains 

0215 

0.2910 

0190 

0  3562 

0  188 

03508 

Canola 

0446 

0.1741 

0.337 

02131 

0  346 

0.2010 

Other  oilseeds 

0  076 

0.6168 

0084 

07551 

0080 

0.7437 

Other  crops 

0.367 

01491 

0.268 

0.1825 

0  269 

0.1797 
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Table  5:    Number  of  statistically  significant  parameters 


Application 


Model  specification 


Number  of  statistically  significant  parameters 


5%  level 


10% 
level 


15% 
level 


20% 
level 


Total 


Application  1 


A1 
B1 
C1 


14 
19 
14 


5 
8 
5 


7 
0 
6 


0 

1 
1 


26 
28 
26 


Application  2 

A2 

17 

6 

7 

0 

30 

B2 

15 

3 

12 

0 

30 

C2 

14 

6 

10 

0 

30 

Note:  The  critical  values  are  1.960,  1.645,  1.439  and  1.281  for  a  5%,  10%,  15%  and  20%  level  of  significance, 
respectively. 


Table  6:    Entropy  ratio  tests  of  model  performance 

Application 

Model  specification 

Entropy  ratio  statistic 

Application  1 

A1 
B1 
C1 

15.284 

7.688 

7.170 

Application  2 

A2 
B2 
C2 

16.624 

9.718 

5.522 

Note:  x"(5)=  11.07  for  a  level  of  significance  of  5%. 


Table  7:    Normalized  entropy  indicators  of  the  ME-tobit  models 

Application 

Model  specification 

S(p) 

S(w) 

Application  1 

A1 
B1 
C1 

0.5463 
0.7375 
0.7706 

0.9033 
0.8983 
0.9003 

Application  2 

A2 
B2 
C2 

0.5110 
0.6999 
0.8363 

0.9529 
0.9469 
0.9458 

4.4     Predictive  power  of  the  ME-tobit  estimates 

The  various  model  specifications  can  also  be  assessed  in  terms  of  their  predictive  power. 
Many  different  indicators  can  be  used  for  this  purpose.  Here,  we  decided  to  use  the 
"pseudo-R2"  which  is  defined  as  the  square  of  the  correlation  coefficient  between  the  pre- 
dicted and  the  observed  values  for  each  cost  item.  This  indicator  is  positive  but  smaller  than 
one.  The  predictive  power  of  the  model  improves  as  the  value  approaches  one. 
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The  "pseudo-R2"  indicator  was  computed  for  each  input,  according  to  the  six  model  specifi- 
cations. A  closer  inspection  of  the  results  presented  in  Table  8  reveals  that  the  values  are 
higher  than  0.800  for  all  the  inputs,  except  two.  For  paid  salaries,  the  "pseudo-R2"  takes  a 
value  of  about  0.501  for  all  model  specifications.  For  fallow  costs,  it  gravitates  around  the 
value  of  0.714. 

Table  8:    Predictive  power  of  the  ME-tobit  models  in  terms  of  the  "pseudo-R2" 
statistic 


Input  cost 

Model  specification 

Application  1 

Application  2 

A1 

B1 

C1 

A2 

B2 

C2 

Seeds 

0.911 

0.908 

0.908 

0.910 

0.907 

0.904 

Fertilizers 

0.906 

0.906 

0.906 

0.894 

0.895 

0.895 

Pesticides 

0.924 

0.920 

0.922 

0.919 

0.915 

0.915 

Other  direct  material  input 

0.854 

0.849 

0.850 

0.854 

0.852 

0.850 

Fuel,  power  and  grease 

0.885 

0.879 

0.882 

0885 

0.882 

0.882 

Repairs 

0.828 

0.818 

0.819 

0.828 

0.818 

0.811 

Paid  salaries 

0.496 

0.497 

0.493 

0.507 

0.504 

0.504 

Other  fixed  cash  expenses 

0.918 

0.914 

0.916 

0.923 

0.920 

0.920 

Fallow  costs 

N/A 

N/A 

N/A 

0.723 

0.714 

0.705 

Net  operating  income 

0.936 

0.932 

0.934 

0.926 

0.924 

0.923 

Note:  "N/A"  means  not  applicable. 

4.5  General  assessment  of  the  ME  estimates 

What  conclusions  can  be  drawn  from  these  estimation  results?  Is  it  justifiable  to  say,  by  look- 
ing at  these  various  indicators  and  statistical  tests,  that  a  particular  model  specification  is 
superior  to  the  others?  If  the  assessment  is  conducted  in  terms  of  statistical  performance  and 
predictive  power,  it  is  rather  difficult  to  rank  the  models  in  terms  of  their  performance.  On 
the  other  hand,  if  one  uses  the  NE  criterion,  more  conclusive  statements  can  be  made.  In  fact, 
one  could  say  without  any  doubt  that  Model  C2,  which  is  characterized  by  an  "asymmetric" 
and  "left-skewed"  spacing  and  a  reduced  interval  of  the  support  values  z,  is  decisively  the 
best  model. 

4.6  Validation  of  the  estimated  cost  allocation  coefficients 

So  far,  the  ME  results  have  been  studied  only  from  econometric  and  statistical  points  of  view, 
without  providing  any  comments  on  (the  plausibility  of)  their  magnitudes.  Evaluating  the 
magnitudes  of  the  coefficients  is  conducted  in  two  stages,  based  on  the  results  of  the  two 
model  specifications  selected,  Models  CI  and  C2.  After  a  first  evaluation  in  general  terms  we 
conducted  a  comparison  of  the  ME  estimates  with  the  available  observed  coefficients  using 
the  cost  data  appearing  in  the  crop  enterprise  database. 

In  general,  most  of  the  ME  estimates  obtained  for  the  direct  inputs,  like  seeds,  fertilizers, 
pesticides,  fuel,  power  and  grease  and  repairs  take  values  smaller  than  0.20.  There  are  a  few 
exceptions  to  this  rule,  however.  For  example,  in  the  case  of  pesticides,  the  coefficient 
associated  with  the  other  oilseeds  output  is  0.386  in  Model  CI  and  0.365  in  Model  C2.  In 
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addition,  a  surprising  result  was  obtained  for  the  coefficients  in  the  equation  for  fertilizers  in 
Application  1  (Models  Al,  Bl  and  CI).  The  Application  1  coefficients  in  the  equations  for 
wheat,  other  grains  and  canola  are  much  higher  than  the  corresponding  coefficients  in  Appli- 
cation 2  (Models  A2,  B2  and  C2).  The  Application  2  coefficient  values,  however,  are  more  in 
line  with  prior  expectations.  Also  the  values  of  the  coefficients  in  the  equation  for  net  operat- 
ing income  changed  drastically  from  Model  CI  to  Model  C2  (the  two  best  model  specifica- 
tions). Specifically,  the  wheat  coefficient  increased  from  0.202  in  Model  CI  to  0.278  in  Model 
C2.  Also,  the  magnitudes  of  the  coefficients  associated  with  other  grains,  canola  and  other 
oilseeds  decreased  significantly  for  these  two  model  specifications.  Overall,  the  ME  estimates 
obtained  with  Model  C2  seem  to  be  more  realistic.  In  addition,  they  are  in  line  with  similar 
estimates  obtained  by  the  authors  of  this  study  for  crops  in  Brittany  (see  Leon  et  al.  1999). 

Furthermore,  it  is  instructive  to  investigate  whether  the  addition  of  fallow  costs  to  the  list  of 
inputs  leads  to  a  change  in  some  coefficients.  A  cross-checking  of  the  coefficients  (reported  in 
Tables  3  and  4)  in  the  equations  for  pesticides,  fuel,  power  and  grease  and  repairs  (inputs 
used  in  fallow-maintenance  operations)  reveals  that  the  coefficients  associated  with  wheat, 
canola  and  other  oilseeds  changed  drastically  from  Model  CI  to  Model  C2.17  A  similar  pat- 
tern was  also  observable  for  the  coefficients  in  the  net  operating  income  equation.  In  most 
other  cases,  the  values  of  coefficients  remained  fairly  stable  or  they  increased  only  margin- 
ally. 

Finally,  the  validity  of  the  ME-tobit  results  is  tested  by  comparing  the  estimated  cost- 
allocation  coefficients  aik  with  the  actual  or  observed  cost  allocations  aik  that  are  available 
from  the  data  set  at  hand  (cost  enterprise  database).  Although  the  actual  cost  allocations  are 
available  at  the  level  of  the  individual  crop  enterprises  (indexed  j  =  1,...,  /),  the  comparison  is 
possible  only  at  the  level  of  the  aggregated  crop  outputs  (indexed  k  =  l,..  .,K),  given  the  nature 
of  the  ME  estimates.  The  comparison  is  carried  out  for  the  five  direct  input  categories  (seeds, 
fertilizers,  pesticides,  fuel,  power  and  grease  and  repairs). 

Before  proceeding  with  this  comparison  however,  weighted  averages  of  the  actual  cost 
allocations  aij{k}  are  also  calculated  (for  information  purposes)  at  the  level  of  the  individual 
crop  enterprises  j{k}  belonging  to  each  aggregated  crop  category  k.  The  weights  are  the 
shares  of  each  individual  farm's  enterprise  acreage  (aCvm)  in  the  corresponding  total 
enterprise  acreage  of  the  sample  \acj(ky  ' 

ac'k)  ac'k) 


and 


SUM,(ac'm)     acm 


<*ij\k)=SUM1{y»)c\J{k))  (25) 


t 


where  cij(k)  is  the  actual  or  observed  cost  for  input  i  allocated  to  crop  enterprise  j  belonging 
to  aggregated  crop  k. 


17.  For  wheat,  the  pesticide  coefficient  decreased  from  0.090  in  Model  CI  to  0.064  in  Model  C2;  for  canola,  the 
corresponding  coefficient  increased  from  0.056  in  Model  CI  to  0.065  in  Model  C2;  for  other  oilseeds,  the 
corresponding  coefficients  decreased  from  0.386  in  Model  CI  to  0.365  in  Model  C2. 
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Next,  weighted  averages  of  the  actual  cost  allocations  at  the  level  of  the  aggregated  crops  k 
(for  k  =  /,...  ,K)  are  calculated,  where  the  weights  are  the  shares  of  each  enterprise's  crop 
acreage  (acjk)  in  the  corresponding  total  crop  acreage  of  the  sample  (ack) : 

acm      SUMt(ac '  k) 

Wjk  =  — —  = J±-L-  (26) 

ack       SUM  (ac  {k)) 


and 


a^SUM^W^a,^)  (27) 


The  calculated,  actual  cost  allocations,  both  at  the  level  of  the  individual  crop  enterprises  and 
the  aggregated  crops  (appearing  under  the  heading  "aggregated  other  crops"),  are  given  in 
Table  9,  along  with  the  corresponding  ME  estimates  at  the  aggregated  level.  Overall,  the  ME 
estimates  for  the  five  direct  inputs  under  consideration  are  clearly  within  the  range  of  the 
observed  coefficients. 

A  closer  inspection  of  the  actual  (calculated)  and  estimated  (ME)  results  in  Table  9  leads  us  to 
make  four  comments: 

•  For  wheat,  which  is  the  main  crop  grown  in  Saskatchewan,  the  average  observed  cost 
allocations  are  smaller  that  the  corresponding  ME  cost-allocation  estimates  for  seeds,  fuel, 
power  and  grease  and  repairs.  For  fertilizers,  the  average  observed  cost  allocation  (which 
is  equal  to  0.175)  is  much  higher  than  the  ME  estimate  generated  by  Model  C2  (0.099),  but 
smaller  than  the  value  obtained  in  Model  CI  (0.193).  For  pesticides,  the  observed  and 
estimated  cost  allocations  are  very  close.  With  regard  to  the  sum,  per  output,  of  the  cost- 
allocation  coefficients  for  the  five  direct  inputs,  the  discrepancy  between  the  observed  and 
estimated  totals  is  smaller  than  10  percent  for  Model  C2. 

•  For  other  grains,  Model  C2  is  reproducing  the  average  observed  cost  allocations 
reasonably  well  for  three  direct  inputs,  namely  seeds,  pesticides  and  repairs.  By  contrast, 
the  cost-allocation  estimates  for  fuel,  power  and  grease  are  completely  out  of  line  with 
their  observed  counterparts.  For  fertilizers,  the  estimate  generated  by  Model  C2  is  very 
close  to  the  observed  cost  allocation.  Finally,  (as  for  wheat)  the  sum,  per  output,  of  the 
estimated  cost-allocation  coefficients  for  the  five  inputs  is  again  very  close  to  the  observed 
total. 

•  The  cost-allocation  estimates  for  canola  are  quite  reasonable  for  three  inputs,  seeds, 
pesticides  and  repairs.  The  cost-allocation  coefficients  in  the  equations  for  the  other  inputs 
were  overestimated  by  a  significant  amount.  Consequently,  it  is  not  surprising  to  find  a 
wide  discrepancy  between  the  sum  of  the  estimated  and  observed  cost  allocations  for  the 
five  direct  inputs. 

•  For  the  remaining  two  crop  categories,  other  oilseeds  and  other  crops,  one  can  discern  a 
wide  discrepancy  between  the  average  observed  and  estimated  coefficients.  This  finding 
is  not  really  surprising,  however,  given  the  heterogeneous  composition  of  these  two  crop 
aggregates. 
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Maximum  entropy  applied  to  the  Saskatchewan  data  sample 


Given  the  empirical  evidence  above,  two  main  findings  emerge: 

•  Using  the  ME  methodology  to  estimate  cost-allocation  coefficients  is  a  valid  approach  in 
the  sense  that  the  estimated  coefficients  fall  within  the  range  of  their  observed 
counterparts. 

•  Some  minor  qualifications  are  needed,  depending  on  the  category  of  crops  under 
consideration. 

Hence,  for  wheat,  other  grains  and  to  a  lesser  extent  also  canola,  the  estimated  and  observed 
cost-allocation  coefficients  are  quite  comparable.  For  the  remaining  crop  categories,  other 
oilseeds  and  other  crops,  some  caution  is  in  order  when  comparing  the  estimated  and 
observed  cost  allocations.  The  discrepancies  found  are  not  really  surprising,  however, 
though  they  may  require  further  investigation  by  creating  more  homogenous  crop  catego- 
ries. 
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Section  5:   Concluding  remarks 

This  study  described  the  results  of  the  application  of  the  Maximum  Entropy  (ME)  methodol- 
ogy to  estimate  the  enterprise-level  cost  allocations,  using  a  sample  of  data  from  30  Saskatch- 
ewan crop  farms.  Although  many  different  methods  could  have  been  used  (and  have  been 
used  by  several  analysts),  this  study  opted  for  the  ME  methodology.  Apart  from  the  fact  that 
the  ME  methodology  offers  a  less  expensive  alternative  to  conducting  a  special  survey  — a 
property  which  the  ME  methodology  has  in  common  with  other  estimation  techniques,  such 
as  OLS  or  LP  — the  choice  for  the  ME  methodology  was  primarily  motivated  by  its  properties 
of  flexibility,  transparency  and  relative  ease  of  implementation.  These  properties  make  it  val- 
uable to  both  practitioners  and  policy  analysts.  Finally,  using  the  ME  methodology  allowed 
us  to  overcome  the  usual  problems  of  implausible  outcomes,  such  as  negative  signs  for  the 
coefficients  (OLS)  and  corner-point  or  zero  solutions  (LP). 

Empirical  results  indicate  that  the  estimated  cost-allocation  coefficients  are  "realistic",  falling 
within  the  range  of  observed  values.  This  finding  holds  true  especially  for  well  defined  crop 
enterprises  (wheat,  other  grains  and  canola)  but  to  a  lesser  extent  for  the  two  heterogeneous 
crop  categories,  other  oilseeds  and  other  crops.  Given  the  exploratory  nature  of  the  present 
study,  further  work  is  needed  to  obtain  more  efficient  and  realistic  estimates  —  for  example, 
by  using  larger  samples,  exploiting  the  "richness"  of  panel  data  and/ or  employing  more  reli- 
able prior  information  on  the  support  values  of  the  cost-allocation  coefficients. 

Some  other  practical  and  methodological  problems  remain.  Specifically,  one  issue  that  needs 
further  investigation  is  how  to  model,  in  a  more  satisfactory  way,  the  decision  process  of  fal- 
low land  maintenance.  This  would  require  the  use  of  sample  data  for  several  years  to  capture 
the  delayed  effect  of  fallow  on  crop  farmers'  returns,  or  the  use  of  more  homogeneous  crop 
farms  —  for  example,  crop  farms  that  are  located  in  the  same  region  or  that  are  using  similar 
technologies.  Another  issue  that  deserves  more  emphasis  has  to  do  with  the  necessity  that 
crop  and  livestock  enterprises  be  described  and  defined  with  great  accuracy. 


To  show  the  merits  of  the  ME  methodology,  we  would  need  to  conduct  a  more  in-depth  comparative 
analysis  of  the  different  estimation  procedures  using  Monte  Carlo  simulation  experiments. 
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Appendix  A:  Definition  of  selected  groups  of 

crops 


Aggregated  crops 


Number 


Name 


Crop 
code 


Enterprise  label 


Product 
code 


Fallow 


1 
2 
201 
17 
52 


Fallow,  Summerfallow(SMF),  Rent  Summerfallow 

Chemical  fallow,  Cons  Fallow 

Green  Manure  sf 

Grass,  IWG  on  rented.  Intermediate  Wheat  GR 

Greenfeed 


1 

1 

1 

16 
35 


Wheat 


103 


10 
110 


43 
143 


Wheat  on  fallow,  Wheat  on  C.F.L,  CRSH-Wheat-Fallow,  CWRS/Fallow 
Wheat/fallow  Rent,  Rent/West,  Rent  East,  Wheat/summerfallow, 
Wheat/stubble,  Wheat/STBL/OWNED,  Canadian  Red  Spring  (CRSH)/ 
wheat-Fall, 

Wheat  on  Wheat,  Wheat  on  canola/stubble,  CSWH-stubble,  CS  durum- 
stubble,  S-wheat/stubble  zero  til,  Rented  wheat/stubble,  CWRS/stubble, 
Tea/Can  St,  Wheat/Rent  stubble,  Wheat/wheat  stubble,  WHT/rent, 
CSWHJ-SJ,  Wheat/stubble/Rent,  Canadian  Prairie  Spring  CPS/fallow, 
CPS  owned,  CPS/stubble, 

Genesis  on  stubble,  Canadian  Prairie  Spnng  (CPS)  on  fallow 

Pr  Spring/on  stubble,  Grandi  on  Canst,  CPS  Wheat/  on  stubble.  CWES/ 
rent/fallow,  CWES/owned/stubble,  CWES/rent/stubble,  CWES/fallow, 
Durum  summer  fallow.  Durum/fallow/rent,  Durum  chemical  fallow  Durum/ 
fallow/rent/FJS,  Durum/fallow(Tref),  Durum  fallow  South 

CRSH  durum/Fallow 

A  Durum/stubble,  Durum/stubble,  Com  Durum  on  Lentils,  Durum/stubble/ 
rent,  Durum  Canadian  Prairie  spring  (CPS)  stubble 


10 
10 


43 
43 


Other  grains 


105 

108 

109 


Barley  on  fallow,  Barley/rented,  Barley/owned,  M/barley/fallow,  Barley/ 
rented/fallow 

Barley/stubble,  Malt  barley  on  stubble,  Feed  barley,  Good  barley, 
Malt  barley/canola  stubble 

Oats,  Oats  on  stubble 

Fall  rye  on  stubble 


IV 


Canola 


106 


Canola  summerfallow.  Rent  canola/fallow,  Canola  on  W  B,  Canola/fallow. 
Canola/fallow/cropshare.  Canola/rented/fallow 

Canola/stubble,  Canola  on  W  B,  Rape  st,  W  Canola  /stubble,  Canola/stub- 
ble owned 


Other 

7 

Flax 

oilseeds 

107 

Flax  on  stubble,  Flax  on  wheat/barley 

111 

Sunola/stubble 

14 

Mustard  fallow,  Y  Mustard  fallow 

114 

Mustard  on  stubble 

7 
11 
14 
14 
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Aggregated  crops 


Number 


Name 


Crop 
code 

Enterprise  label 

Product 
code 

153 

Certified  wheat  fallow,  Certified  wheat/chemical  fallow,  Pedigree  wheat  on 
fallow 

3 

157 

Pedigree  flax  summerfallow 

11 

12 

Canary  seeds,  Canary  seeds/stubble,  Canary  seeds/fallow 

12 

112 

Canary  seeds/stubble 

12 

13 

Lentils,  Lentils/fallow 

13 

113 

Lentils  (Rape  St),  Lentils/canola  stubble,  Lentils/stubble/rented,  Lentils  on 
wheat  stubble,  CRSH  Lentils-St,  Lentils/stubble,  Richlea/stubble,  Richlea 
lentils,  Lentils  rent,  Lentils/stubble  Owned,  Lentils/Canola  stubble,  Cert  Est 
Lentils/fallow 

13 

163 

Pedigree  Lentils/wheat  stubble 

13 

183 

Peas 

13 

15 

Peas/stubble,  peas/cereal  stubble,  Peas/stubble/rent,  Peas  on  W  B. 

15 

115 

Crsh  Peas/Stubble 

15 

17 

Dahurian 

16 

17(1) 

Alfalfa(rent) 

16 

19 

Alfalfa  hay 

16 

20 

Seed  red  clover 

23 

16 

Hay 

33 

50 

Pasture 

33 

53 

Chickpeas/owned/fallow,  Triticale,  Caraway,  Coriander/fallow,  Coriander/ 
stubble, 

36 

45 

Borage/stubble 

45 

193 

Certified  Durum,  Fallow 

43 

213 

Ped  Durum/dir  stubble 

43 

VI 


Other  crops 
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Appendix  B:  Data  set 


Table  B-1:  Data  set  -  Net  returns 


Farm 

Net  returns 

Identification  No. 

Wheat 

Other  grains 

Canola 

Other  oilseeds 

Other  crops 

Y1 

Y2 

Y3 

Y4 

Y5 

25 

13320.00 

31065  00 

5381685 

5600  00 

69766  80 

27 

387450.00 

3196800 

372937.50 

132600.00 

100270.5 

30 

89082.75 

44200  00 

91650  00 

0.00 

161  00 

33 

3661.51 

12576.61 

33690.10 

0.00 

22180  500 

35 

129365  60 

86269  20 

210754.20 

0  00 

80821.40 

42 

18620  93 

20550.26 

9253585 

000 

32980.57 

56 

218625  00 

131800  00 

287087  50 

29448  00 

05320  00 

58 

65827  26 

0  00 

0.00 

0.00 

75450.37 

77 

31500  00 

28644  00 

000 

11700  00 

000 

98 

40262  30 

0.00 

0.00 

34220  00 

53722.50 

99 

5957360 

6864.25 

4156.56 

0.00 

2395.98 

109 

72168  00 

0.00 

0.00 

0.00 

367590  00 

111 

156172  40 

9584.50 

000 

7424.00 

7739200 

121 

66494.00 

18320.00 

46800  00 

000 

0.00 

123 

12277400 

0.00 

12990682 

0.00 

179622.00 

129 

10991.59 

20784  09 

23120.76 

0  00 

11073  04 

136 

174489  90 

0.00 

312841  32 

4327960 

138892  600 

167 

163335  68 

8178.19 

000 

0.00 

17139  60 

182 

102376  56 

437500 

50920  00 

0.00 

43417.80 

185 

176017  50 

21562.50 

206009  78 

0  00 

81366.00 

186 

57750  00 

0.00 

0.00 

0.00 

15120.00 

201 

79780  00 

110656  00 

000 

0.00 

64631  13 

215 

48737.70 

0  00 

80850  00 

15600  00 

1417500 

218 

9508924 

621000 

82881  25 

000 

49600  00 

220 

000 

32152  50 

47190  00 

6300  00 

1381250 

221 

20312.50 

107306  50 

35640.00 

0.00 

0.00 

259 

30624  00 

52195.50 

65483.31 

1817200 

000 

286 

18564.26 

0.00 

7933998 

0  00 

19264  00 

289 

76238  60 

0.00 

56880  30 

904500 

0  00 

295 

48915.90 

1281800 

24360  00 

16652  34 

0  00 

Farm  cost  allocation  based  on  the  Maximum  Entropy  methodology 


B-37 


Appendix  B 


Table  B-2:  Data  Set  -  Costs 


Costs 

Farm 
Identification 

Other 

direct 

material 

inputs 

Fuel, 

Paid 
salaries 

Otherfixed 

Net 

No. 

Seeds 

Fertilizers 

Pesticides 

power 

and 

grease 

Repairs 

cash 
expenses 

operating 
income 

X1 

X2 

X3 

X4 

X5 

X6 

X7 

X8 

X9 

25 

1294045 

20784  00 

17932  90 

6717  00 

14298  00 

19086  00 

3000  00 

4219900 

36611  30 

27 

59645.38 

95770.35 

145437  220 

32309  00 

6267000 

6620000 

36500  00 

172018.00 

354676.05 

30 

10805  00 

2224269 

16957.76 

5738.00 

21943.00 

18762.00 

7000.00 

40048.00 

81597  31 

33 

6971.93 

10205.60 

8999  80 

5159.00 

7862.88 

6321.24 

0.00 

23688.00 

2900.27 

35 

2337436 

32804.00 

26712.50 

3095900 

20754.00 

19577.00 

138264  00 

80341.00 

134424.55 

42 

13683.99 

24425.75 

26505  00 

6984.00 

11772.00 

13244.00 

400000 

32221.00 

31851  87 

56 

1759925 

50617  005 

28246  20 

1491800 

4099500 

2807600 

55670.00 

12964500 

40651406 

58 

1530  00 

0.00 

17437.80 

5154.00 

3185.19 

13571.00 

1428.00 

23120.00 

75851.64 

77 

6545.50 

11400  00 

19636  00 

3970.00 

7904.25 

7453.00 

0.00 

10452.00 

4483.25 

98 

18145.51 

367  20 

21727.53 

7505.50 

21475.00 

6186.00 

12500.00 

52329.00 

-1203093 

99 

981.00 

4688.33 

8753  31 

5952.00 

472065 

6938  00 

668200 

34734  00 

-458.90 

109 

23953.75 

18826.00 

43070  58 

22704  50 

11536  00 

17646  00 

12400.00 

38288.00 

251333.18 

111 

20933.00 

2146065 

27966  00 

11763.00 

12869  00 

3065200 

5640.00 

64692.00 

54597  25 

121 

6720.08 

12612.20 

15822.45 

13429  00 

14055.00 

15394  00 

3500.00 

49697  00 

384.27 

123 

28519.55 

34547.04 

26045.07 

25870.00 

12780.00 

8150  00 

15000.00 

51756  00 

229635  16 

129 

7888  75 

2721.95 

4324.50 

15029  00 

0  00 

0.01 

5400.00 

23920  00 

6685.27 

136 

34526.91 

83023  10 

91816.31 

21264.00 

24306  00 

10075.00 

16500  00 

69123.00 

318869  10 

167 

13187.49 

10809  00 

22969.18 

10362  00 

17658  00 

1410600 

38386  00 

3881000 

22365.79 

182 

17880  04 

11079.12 

12208  98 

7701.00 

14577.00 

26467.00 

33500  00 

20155.00 

57521.22 

185 

39402.88 

28731  00 

4510850 

20554.00 

22920.00 

19500.00 

40000  00 

3871800 

230021.41 

186 

5245.00 

9807.00 

6555.40 

3900.00 

1368.00 

55501 

0.00 

2510.00 

42929.59 

201 

10337.60 

5674.50 

7568  00 

7406.00 

3911.66 

3531  00 

380  00 

23018.00 

189240  37 

215 

888671 

9926.25 

9608.21 

8996.50 

12997  00 

5547.00 

12000.00 

29497.00 

61904.04 

218 

18406.25 

36251  10 

11112.05 

7268.00 

13692.00 

20344  00 

7000  00 

25752.00 

93955.09 

220 

6216.25 

16857.50 

15980.13 

6516.00 

1301400 

11431.00 

7000.00 

26062.00 

-3621.88 

221 

499388 

23817  00 

7624.20 

9871.00 

13814.00 

11000  00 

4500.00 

31667.00 

55971.93 

259 

6350.00 

18384.00 

15793.20 

4757.00 

8567.35 

9700.00 

2511.00 

20483.00 

79929.26 

286 

7281.00 

5932.00 

9130.00 

8476.00 

4288.00 

6623.00 

0.00 

18938.00 

56500.24 

289 

3744  00 

2721750 

7455.75 

14845.00 

11002.00 

14681.00 

0.00 

28351  00 

34867.65 

295 

5115.25 

14851.5 

118760 

1450.00 

6571.25 

3910.00 

0.00 

18796.00 

40176.19 
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Appendix  C:  GAMS  program  to  estimate  the 

cost-allocation  coefficients 


*  Filename  "saskat.gms" 


*  The  objective  of  this  GAMS  program  is  to  estimate  cost-allocation  coefficients 

*  for  Saskatchewan  crop  farms.  The  program  develops  the  ME-tobit  specifications 

*  without  fallow  cost  and  with  11  support  values  symmetrically  spaced  over  the 

*  [0,  1]  interval. 


SETS 

L 

raw  index 

T 

farm  index 

I 

input  index 

/  1  *  15  / 

/  1  *  30  / 

/  1  seeds  cost 

2  fertilizers  cost 

3  pesticides  cost 

4  other  direct  material  cost 

5  fuel,  power  and  grease  cost 

6  repairs  cost 

7  paid  salaries  cost 

8  other  fixed  cash  expenses  cost 

9  net  operating  income  / 


K     output  index  /  1  wheat  output 

2  other  grains  output 

3  canola  output 

4  other  oilseeds  output 

5  other  crops  output  / 


M 


index  parameter  support  /  1  *  11  / 


N 


index  error  support  7  1*3/ 
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acres  index   /  1  fallow  land 

2  wheat 

3  other  grains 

4  canola 

5  other  oilseeds 

6  other  crops  / 


$INCLUDE  "saskatl.prn"; 
$INCLUDE  "acres. prn" 

ALIAS  (I,  J)  ; 
ALIAS  (T,S)  ; 

*  Display  the  raw  data  on  input  used,  output  used  (Table  MN) 

*  and  acreages  (Table  AA) 

DISPLAY  MM; 
DISPLAY  AA; 

VARIABLES 
OBJ; 

PARAMETERS 

*  Z (M)  is  the  vector  of  support  values  for  the  input/output  coefficients 


Z  (M)    / 

1 

0 

0 

2 

0 

1 

3 

0 

2 

4 

0 

3 

5 

0 

4 

6 

0 

5 

7 

0 

6 

8 

0 

7 

9 

0 

8 

10 

0 

9 

11 

1 

0/ 

X{T,I) 

XYZ (T, I) 

Y(T,K) 

XR(T) 

FID(T) 

DIF(T) 

DIFF(T) 

V(N,I) 

SIGMA(I, 

J) 

RSQUARE(I) 

RSQUARR(I) 

CRITP 

CRITU 
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GAMS  program  to  estimate  the  cost-allocation  coefficients 


XHAT(T,I) 
XF(T, I) 
W(T,I,N) 
AHAT(K, I) 
AHATU(K,I) 
SUMAHAT (K) 
SUMAHATU(K) 


*   Defining  input  and  output  variables 


FID(T) 

Y(T, "1" 

Y(T, "2" 

Y(T,"3" 

Y(T,  "4" 

Y(T,"5" 

X(T, "1" 

X(T,"2" 

X(T,"3" 

X(T,"4" 

X(T,"5" 

X(T,"6" 

X(T,"7" 

X(T, "8" 

X(T,"9" 


MM(T, 
MM(T, 
MM(T, 
MM(T, 
MM(T, 
MM(T, 
MM(T, 
MM(T, 
MM(T, 
MM(T, 
MM(T, 
MM(T, 
MM(T, 
MM(T, 
MM(T, 


"1") 

"2") 

"3") 

"4") 

"5") 

"6") 

"7") 

"8") 

"  9  " ) 

"10") 

"11") 

"12") 

"13") 

"14") 

"15") 


XR(T)=SUM(K, Y(T,K)  ) -X (T, "1" ) -X (T, "2" ) -X (T, "3" ) -X (T, "4 " ) -X (T, "5" ) -X (T, "6" ) -X  (T, "7' 

X(T,"8")  ; 

DIF(T)=XR(T)-X(T,  "9")  ; 

DISPLAY  FID  ; 

DISPLAY  DIF; 

X  ( T , " 9 " ) =XR ( T ) ; 

*  Estimate  of  the  standard  error  of  the  residuals 

TABLE  VI (N, I) 

12345                      67                        8  9 

1      -6087    -28599.7    -10036.9    -5620.5    -18715.1    -8056.2    -49891.7    -15351.7  -156467.3 

200                      000000  0 

3         6087    -28599.7      10036.9      5620.5       18715.1      8056.2      49891.7       15351.7  156467.3 


*   Establishing   the   "three-sigma"    rule   for   the   residuals 

V(N, I)=3*V1 (N,  I) ; 
DISPLAY   V; 
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POSITIVE  VARIABLES 

*  Probabilities  associated  with  support  values  of  the 

*  cost-allocation  coefficients  and  residuals 

P(K,I,M) 
W1(T,I,N) 

W2(T,I,N) ; 

*  Defining  variables  that  are  positive  or  equal  to  zero 

XYZ(T,I)  =  0$(X(T,I)  LE  0)  +  X (T, I ) $ (X (T, I )  GT  0)  ; 

EQUATIONS 
OBJECTIVE 
ADD1 (K,I) 
ADD21 (T,  I) 
ADD22 (T,I,N) 
ADD31 (T, I) 
ADD32 (T,I,N) 
CON1 (T,I) 
CON2 (T, I) 
SUMRESTR(K) 

*  Objective  function  of  the  ME  problem 

OBJECTIVE..  OBJ  =E=-SUM(I,  SUM(K,  SUM(M,  P (K, I ,M) *LOG ( 1 . e-4+P (K, I , M) ) ) 

-SUM (I,  SUM(T$ (XYZ (T, I)  GT  0), 

SUM(N,  Wl (T,I,N) *L0G(l.e-4+Wl (T, I,N) ) ) ) ) 
-SUM (I,  SUM(T$ (XYZ (T,I)  LE  0), 

SUM(N,  W2 (T,I,N)*LOG(l.e-4+W2 (T,I,N) ) ) ) ) ; 

*  SET  OF  EQUATIONS 

*  Conditions  that  the  sum  of  probabilities  P(.),  Wl ( . )  and  W(.) 
is  equal  to  one 

ADD1(K, I)..    SUM(M,  P(K,I,M))  =E=  1.0  ; 

ADD21 (T,I)$ (XYZ (T, I)  GT  0)..   SUM(N,  W1(T,I,N))  =E=  1.0  ; 
ADD22 (T,I,N)$ (XYZ (T, I)  LE  0)..   W1(T,I,N)  =E=  0.0  ; 

ADD31 (T,I) $ (XYZ(T, I)  LE  0)..   SUM(N,  W2(T,I,N))  =E=  1.0  ; 
ADD32 (T, I,N) $ (XYZ (T,I)  GT  0)..   W2(T,I,N)  =E=  0.0  ; 

*  Input  demand  equations 

CON1 (T,I)$ (XYZ(T, I)  GT  0) . .   XYZ(T,I)  =E=SUM(K,  SUM(M, 

P(K,I,M) *Z(M) ) *Y(T,K) ) 

+SUM(N,  Wl (T,I,N)*V(N,I) )  ; 
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C0N2 (T, I) $ (XYZ(T,I)  LE  0)..   SUM(K,  SUM (M,  P (K, I,M) *Z (M) ) *Y (T, K) ) 

+  SUM(N,  W2 (T, I,N) *V(N, I) )  =L=  0  ; 

*  Adding-up  condition  on  the  cost-allocation  coefficients 
SUMRESTR(K) . .   SUM(I,  SUM (M,  P (K, I ,M) *Z (M) ) )  =E=  1  ; 

*  Resolution  of  program  using  MINOS 

OPTIONS  ITERLIM  =  20000  ; 

OPTIONS  RESLIM=3000; 

MODEL  BRET  /ALL/; 

OPTION  NLP=MINOS5; 

SOLVE  BRET  USING  NLP  MAXIMIZING  OBJ  ; 

*  Display  results 

*  Cost-allocation  coefficients 

AHAT(K, I)  =  SUM(M,  P. L (K, I ,M) *Z (M) )  ; 
DISPLAY  AHAT  ; 

SUMAHAT(K)  =  SUM ( I ,  AHAT (K, I))  ; 
DISPLAY  SUMAHAT  ; 

*  Normalized  entropy  ratios 

CRITP  =  -SUM (I,  SUM(K,  SUM(M,  P. L (K, I ,M) *LOG ( 1 . e-4+P. L (K, I ,M) ) ) ) ) ; 
CRITU  =  -SUM (I,  SUM(T$ (XYZ(T,I)  GT  0), 

SUM(N,  W1.L(T,I,N) *LOG(l.e-4+Wl.L(T,I,N) ) ) ) ) 
-SUM (I,  SUM(T$ (XYZ (T,I)  LE  0), 

SUM(N,  W2.L(T,I,N)+LOG(l.e-4+W2.L(T,I,N) ) ) ) ) ; 
DISPLAY  CRITP; 
DISPLAY  CRITU; 

*  Estimated  residuals 

W(T,I,N)=W1.L(T,I,N)$ (XYZ (T,I)  GT  0)  +  W2 . L (T, I, N) $ (XYZ (T, I )  LE  0); 

*  Estimated  variance/covariance  of  residuals 

SIGMA (I, J)  =  SUM(T, (SUM(N,  W (T, I , N) *V (N, I ) ) 

*  SUM(N,  W(T, J,N)*V(N, J) ) ) )/ (CARD (T) -5) ; 

DISPLAY  SIGMA; 

*  Predicted  values  of  input  demand 

XHAT(T,I)$ (XYZ(T,I)  GT  0)  =  SUM(K,  AHAT (K, I ) * Y (T, K) $ (XYZ (T, I )  GT  >0)); 
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*  Compute   R-square 

RSQUARR(I)     =    SUM(T$(XYZ(T,I)     GT    0),    XHAT (T, I ) *X (T, I ) ) 

*    SUM(T$ (XYZ (T,I)     GT    0),    XHAT (T, I ) *X (T, I ) ) / 

(SUM(T$ (XYZ(T,I)  GT  0),  XHAT (T, I ) *XHAT (T, I ) ) 
*  SUM(T$(XYZ(T,I)  GT  0),  X (T, I ) *X (T, I ) ) ) ; 
DISPLAY  RSQUARR; 

*  Save  estimated  matrix  of  variance-covariance  of  the  estimated 

*  residuals 

FILE  SASKASOL  / "saska . sol"/ ; 
PUT  SASKASOL; 
SASKASOL. PC=5; 
SASKASOL. PW=25  0; 
LOOP ( I , 
LOOP (J,  PUT   SIGMA (I, J)  ; 

); 

PUT  /; 
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variance-covariance  matrix  of  the 
ME -estimated  cost-allocation 
coefficients 


OPTIONS  MEM0RY=12; 

SET  N=30; 

SMPL  1,N; 

LOAD(file='SASKATl.prn' ) FARMID  Yl  Y2  Y3  Y4  Y5  XI  X2  X3  X4  X5  X6  X7  X8  X9; 

GENR  Y9R=Y1+Y2+Y3+Y4+Y5-X1-X2-X3-X4-X5-X6-X7-X8; 

■)    ************************* 
?  DEFINITION  OF  VARIABLES 

9   ************************** 

7 

?  Yl=  WHEAT  OUTPUT 

?  Y2=  OTHER  GRAINS  OUTPUT 

?  Y3=  CANOLA  OUTPUT 

?  Y4=  OTHER  OILSEEDS  OUTPUT 

?  Y5=  OTHER  CROPS  OUTPUT 

7 

?  Xl=  SEEDS  COST 

?  X2=  FERTILIZERS  COST 

?  X3=  PESTICIDES  COST 

?  X4=  OTHER  DIRECT  MATERIAL  INPUTS  COSTS 

?  X5=  FUEL,  POWER  AND  GREASE  COSTS 

?  X6=  REPAIRS  COSTS 

?  X7=  PAID  SALARIES  COSTS 

?  X8=  OTHER  FIXED  CASH  EXPENSES  COSTS 

?  X9R=  NET  OPERATING  INCOME  I 
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SMPL  1,4  5; 

LOAD  RR1;  1000010000100001000010000 
10000100001000010000; 

LOAD  RR2;  0100001000010000100001000 
01000010000100001000; 

LOAD  RR3;  0010000100001000010000100 
00100001000010000100; 


LOAD  RR4;  0001000010000100001000010 
00010000100001000010; 

LOAD   RR5;  0000100001000010000100001 
00001000010000100001; 

MMAKE  RT  RR1  RR2  RR3  RR4  RR5 ; 
MAT  R=RT ' ; 

SMPL  1  30; 

MMAKE  YY  Yl  Y2  Y3  Y4  Y5; 

MAT  YYT=YY ' ; 

MAT  YYYY=YYT*YY; 

SMPL  1  9; 

?  SIGMA1   variance/covariance  matrix  of  the  estimated  residuals 

READ (FILE=' sigmal.prn' ) SIGH  SIGI2  SIGI3  SIGI4  SIGI5  SIGI6  SIGI7  SIGI8  SIGI9; 

MMAKE  SIGMA  SIGH  SIGI2  SIGI3  SIGI4  SIGI5  SIGI6  SIGI7  SIGI8  SIGI9; 

MAT  SIGINV=SIGMA"; 

MAT  II=IDENT(30) ; 

MAT  DDD=SIGINV#YYYY; 

MAT  OMEGA=DDD"; 

MAT  NNINV=(R*OMEGA*RT) " ; 

MAT  OMEGAME=OMEGA-OMEGA*RT*NNINV*R*OMEGA; 

MAT  VAR=DIAG(OMEGAME) ; 

MAT  SE=VAR**0.5; 

PRINT  SE; 

STOP; 

END; 
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Appendix  E:  Normalized  entropy  indicators 

associated  with  the  estimated 
cost-allocation  coefficients 


Model  A1 

Model  B1 

Model  C1 

Model  A2 

Model  B2 

Model  C2 

1.  Seeds 

Wheat 

0470 

0642 

0  704 

0459 

0635 

0813 

Other  grains 

0  171 

0  348 

0  428 

0  189 

0  367 

0606 

Canola 

0317 

0476 

0  546 

0.319 

0480 

0674 

Other  oilseeds 

0476 

0699 

0.738 

0476 

0689 

0865 

Other  crops 

0463 

0672 

0.719 

0455 

0659 

0  844 

2.  Fertilizers 

Wheat 

0  780 

0971 

0  967 

0535 

0782 

0937 

Other  grains 

0716 

0.918 

0  927 

0644 

0872 

0979 

Canola 

0.797 

0.983 

0.976 

0.648 

0.895 

0.984 

Other  oilseeds 

0548 

0771 

0  797 

0477 

0698 

0872 

Other  crops 

0612 

0839 

0  857 

0584 

0835 

0.959 

3.  Pesticides 

Wheat 

0  542 

0  794 

0.797 

0477 

0.677 

0838 

Other  grains 

0  086 

0  154 

0  307 

0  044 

0.080 

0265 

Canola 

0429 

0597 

0  669 

0451 

0654 

0843 

Other  oilseeds 

0979 

0  901 

0.954 

0.972 

0.909 

0730 

Other  crops 

0  567 

0762 

0810 

0.588 

0807 

0  952 

4.  Other  direct  material  inputs 

Wheat 

0318 

0  460 

0.541 

0327 

0478 

0670 

Other  grains 

0382 

0  584 

0643 

0  376 

0570 

0775 

Canola 

0.285 

0  430 

0  503 

0283 

0431 

0616 

Other  oilseeds 

0303 

0522 

0.587 

0294 

0  504 

0  736 

Other  crops 

0  447 

0647 

0  698 

0444 

0641 

0826 

5.  Fuel,  power  and 

grease 

Wheat 

0  421 

0617 

0685 

0440 

0  652 

0  849 

Other  grains 

0617 

0857 

0  862 

0594 

0  825 

0  955 

Canola 

0448 

0679 

0713 

0440 

0677 

0  857 

Other  oilseeds 

0482 

0706 

0  742 

0472 

0685 

0  863 

Other  crops 

0451 

0  701 

0  729 

0439 

0  677 

0  863 
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Model  A1 

Model  B1 

Model  C1 

Model  A2 

Model  B2 

Model  C2 

6.  Repairs 

Wheat 

0545 

0.720 

0.769 

0  546 

0729 

0881 

Other  grains 

0483 

0721 

0.754 

0.459 

0.686 

0871 

Canola 

0202 

0  378 

0.454 

0.188 

0.361 

0.587 

Other  oilseeds 

0.502 

0  720 

0753 

0  502 

0.706 

0872 

Other  crops 

0314 

0.517 

0579 

0303 

0495 

0715 

7.  Paid  salaries 

Wheat 

0504 

0761 

0790 

0548 

0821 

0.956 

Other  grains 

0.738 

0943 

0  940 

0.725 

0.925 

0.996 

Canola 

0621 

0899 

0878 

0  628 

0912 

0.989 

Other  oilseeds 

0458 

0697 

0733 

0.446 

0674 

0859 

Other  crops 

0668 

0925 

0909 

0649 

0900 

0  986 

8.  Other  fixed  cash  expenses 

Wheat 

0785 

0958 

0965 

0779 

0  960 

0999 

Other  grains 

0807 

0978 

0.975 

0751 

0943 

0999 

Canola 

0.642 

0894 

0890 

0621 

0872 

0978 

Other  oilseeds 

0632 

0833 

0847 

0.624 

0815 

0.937 

Other  crops 

0.580 

0833 

0841 

0559 

0  800 

0.942 

9.  Fallow  cost 

Wheat 

N/A 

N/A 

N/A 

0249 

0  402 

0.609 

Other  grains 

N/A 

N/A 

N/A 

0  565 

0783 

0.923 

Canola 

N/A 

N/A 

N/A 

0209 

0.371 

0.574 

Other  oilseeds 

N/A 

N/A 

N/A 

0  392 

0  594 

0.793 

Other  crops 

N/A 

N/A 

N/A 

0255 

0440 

0.659 

10.  Net  operating  income 

Wheat 

0.788 

0982 

0.974 

0902 

0.995 

0918 

Other  grains 

0.807 

0.977 

0973 

0.802 

0969 

0  998 

Canola 

0.955 

0978 

0.988 

0  986 

0934 

0782 

Other  oilseeds 

0506 

0.751 

0.770 

0499 

0736 

0.892 

Other  crops 

0.939 

0995 

0997 

0940 

0997 

0.930 
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